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SUMMARY 


The  prototype  instrument  to  determine  the  orientation  of 
a  projectile  in  free  flight  (i.  e.  roll,  pitch  and  yaw 
angles)  was  developed  at  Boeing  and  used  to  perform  tests  at 
the  U.S.  Army  Ballistic  Research  Laboratory  (BRL)  at 
Aberdeen  Proving  Ground,  MD.  The  instrument  incorporated  two 
visible  gas  lasers  (to  provide  a  two  color  collimated 
illuminating  beam)  and  a  pair  of  galvanometer  beam  deflectors 
to  generate  a  circular  angular  dither  of  the  output  beam. 
All  three  angles  can  be  determined  from  the  detected  light 
reflected  by  a  retroref lector/hologram  combination  mounted  on 
the  projectile.  The  prototype  instrument  was  designed  to 
make  three-axis  angle  measurements  on  in-flight  or  in-bore 
spinning  or  non-spinning  projectiles,  and  to  make  two-axis 
angle  measurements  of  gun  tube  motion. 

Tests  were  performed  by  both  Boeing  and  BRL  personnel  at 
the  Aerodynamics  Range  at  BRL  during  a  two  week  period  in 
August,  1985.  The  test  instrumentation  was  then  left  at  BRL 
for  additional  testing  by  BRL  personnel.  A  summary  of  the 
test  results  follows:  1)  All  three  angles  of  a  spinning 
projectile  can  be  measured  with  a  fixed  beam  and  two 
detectors,  where  one  detector  senses  fringe  amplitude  and  the 
other  detector  senses  roll,  2)  The  dither  rate  (3  kilohertz) 
was  too  low  for  the  non-spinning  projectiles  tests,  3)  No 
in-bore  data  was  obtained  due  to  obturation,  4)  Dynamic 
measurements  of  both  angle  and  position  of  a  gun  tube  were 
demonstrated . 

As  a  result  of  the  test,  the  following  problem  areas 
have  been  identified:  1)  The  test  results  show  that  a  system 
with  a  higher  dither  frequency  must  be  developed,  2)  The  use 
of  two  independent  lasers  complicate  the  alignment  procedure. 
It  would  be  more  reasonable  to  use  a  two  color  laser  (e.g. 
argon  ion  laser),  3)  Some  of  the  data  contains  an  intensity 
variation  at  the  fundamental  of  the  dither  frequency.  This 
could  be  due  to  the  projectile  lying  at  the  edge  of  the 
"overlap''  of  the  dithered  beams.  A  larger  beam  should 
eliminate  this  problem  and  simplify  the  alignment  procedure, 
4)  The  reflection  coefficient  versus  angle  was  assumed  to  be 
sinusoidal  for  the  present  data  reduction  algorithm.  In 
practice,  the  reflection  coefficient  is  not  sinusoidal,  and 
the  shape  of  the  function  changes  with  relative  humidity  and 
laser  wavelength.  Additional  development  is  needed  on  an 
algorithm  based  on  a  more  realistic  analytic  model. 


INTRODUCTION 


This  document  is  the  final  report  on  Phase  1  of  contract 
DAAK1 1-84-C-0095 .  The  contract  is  the  first  part  of  a 
program  whose  goal  is  to  develop  hardware  for  measuring  the 
angular  orientation  of  projectiles  in-flight  and  in-bore, 
and  the  angular  response  of  gun  tubes  during  firing.  The 
technology,  originally  developed  for  measuring  the  pitch 
angle  of  models  in  the  Boeing  Transonic  Wind  Tunnel,  was 
shown  to  be  applicable  to  spinning  projectile  yaw  angle 
measurement  in  a  prior  test  program  (Reference  2)  at  the  U.S. 
Army  Aerodynamics  Range.  Significant  improvements  to  the 
measurement  techniques  were  made  in  the  present  contract. 

A  two-laser  breadboard  system  was  assembled  at  Boeing 
for  evaluation  at  the  Ballistic  Research  Laboratory.  The 
breadboard  was  designed  to  make  3-axis  angle  measurements  on 
in-flight  and  in-bore  non-spinning  or  spin-stabilized 
projectiles,  and  2-axis  angle  measurements  on  gun  tubes.  For 
spin-stabilized  projectiles,  three-axis  angle  measurements 
were  demonstrated.  Signals  were  obtained  from  non-spinning 
projectiles,  but  could  not  be  converted  into  meaningful 
information  due  to  the  inability  of  the  breadboard  optical 
system  to  measure  the  orientation  of  projectiles  possessing 
high  yawing  rates.  No  in-bore  data  was  obtained  due  to 
degradation  of  the  incident  laser  beam  quality  resulting  from 
the  leakage  of  hot  combustion  gases  around  the  projectile. 
Dynamic  measurement  of  both  angle  and  position  were 
demonstrated  for  the  gun  tube  tests. 

The  report  contains  a  description  of  the  breadboard 
system,  the  modifications  required  for  the  various 
measurement  tasks,  the  test  conditions,  and  the  test  results. 
Data  reduction  algorithms  developed  by  Boeing  are  in  Appendix 
B.  Preliminary  design  information  for  Phase  II  hardware, 
based  on  the  results  of  this  evaluation,  are  included. 
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The  test  setup  consists  of  an  optical  unit  and 
supporting  electronics  located  in  the  range  together  with 
data  acquisition  and  processing  equipment  placed  in  the 
control  room.  The  optics  will  be  discussed  first. 

Figure  1  shows  the  arrangement  of  optical  elements  on  a 
standard  0.6x1. 2  meter  setup  board.  Red  and  blue  lasers  are 
used  to  provide  orthogonal  angle  measurements  on  the 
projectile.  The  beams  are  combined  in  a  dichroic  cube 
beamsplitter  and  then  directed  to  the  filter-splitter.  A 
lens  in  the  red  beam  is  used  to  adjust  the  beam  diameter  at 
the  filter-splitter  so  that  the  red  and  blue  beams  have  the 
same  diameter  when  they  exit  the  unit. 

The  filter-splitter  is  a  Boeing  device  (U.  S.  Patent  No. 
4,329,059)  that  spatially  filters  the  outgoing  beam  while 
transmitting  most  of  the  return  beam  energy  to  the 
photodetectors .  The  coaxial  red  and  blue  beams  at  the  input 
to  the  filter-splitter  are  focused  by  a  microscope  objective 
lens  onto  a  small  elliptical  mirror  on  a  tilted  glass  plate 
which  spatially  filters  the  reflected  beams  and  insures  that 
the  red  and  blue  beams  are  collinear.  The  elliptical  mirror 
is  small  compared  to  the  diffraction  pattern  of  the  return 
beams  from  the  reflector  on  the  projectile,  so  more  than  90% 
of  the  return  energy  passes  on  to  the  photodetectors. 

A  two-axis  galvanometer  type  beam  deflector  generates  a 
circular  dither  (or  conical  scan).  The  relay  lens  forms  an 
image  of  the  filter-combiner  mirror  in  the  focal  plane  of  the 
large  collimating  lens.  The  function  of  the  relay  lens  is  to 
reduce  lateral  motion  of  the  spot  at  the  projectile.  (The 
translation  is  zero  in  the  plane  where  the  combination  of  the 
collimating  and  relay  lenses  forms  an  image  of  a  galvanometer 
mirror. )  The  output  beam  is  collimated  and  reflected  into 
the  projectile  path  with  a  pair  of  mirrors  (not  shown). 

The  reflector  assembly  consists  of  a  hologram  and  a  cube 
corner  reflector,  as  described  in  Reference  2,  with 
modifications  that  allow  3-axis  angle  measurement.  The 
hologram  has  a  central  circular  holographic  grating 
surrounded  by  an  annular  grating  with  equal  area  but 
orthogonal  fringe  orientation.  Red  and  blue  spectral  filters 
that  match  the  hologram  profiles  are  sandwiched  between  the 
hologram  and  the  reflector. 
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DATA  ACQUISITION  ELECTRONICS 


Figure  2  is  a  block  diagram  of  the  data  Acquisition 
Electronics  (DAE).  The  DAE  was  used  to  collect  range  or 
calibration  data.  The  photodetectors  were  Hamamatsu  R761 
photomultipliers  (PMT),  a  head-on  type  with  a  half  inch 
(13mm)  diameter  photocathode.  The  PMT  has  ten  dynode  stages, 
typical  current  gain  of  5.8X10E+5,  and  a  spectral  response 
from  185nm  to  850nm,  peaking  at  420nm.  The  light  entering 
the  PMT’s  was  optically  filtered  for  the  two  optical 
wavelengths  before  application  to  the  PMT's. 

Both  the  red  and  blue  signal  channels  are  identical  with 
the  exception  of  the  overall  gain.  Load  resistors  for  the 
PMT’s  were  10k  ohms.  At  the  optical  assembly  each  PMT  output 
was  amplified  by  a  HP465A  amplifier  with  a  gain  of  20db  and 
filtered  by  a  Krohn-Hite  3342  filter/amplifier.  One  section 
of  the  Krohn-Hite  filter  is  used  as  a  low  pass  filter  with  a 
cutoff  frequency  of  100Hz  and  a  gain  of  20db.  The  low  pass 
section  drives  the  high  pass  section  resulting  in  a  bandpass 
with  cutoff  frequencies  of  100Hz  and  24KHz.  The  high  pass 
section  of  the  Krohn-hite  filter  drives  approximately  100 
feet  of  terminated  coax  cable  between  the  optical  assembly 
and  the  control  room.  The  variable  gain  capability  of  the 
HP467A  was  used  to  adjust  the  input  amplitude  to  the  tape 
recorder  and  analog-to-digital  converter  to  approximately  one 
volt  RMS. 

All  timing  for  the  system  is  synchronously  derived  from 
a  very  stable  frequency  source,  an  HP8460B  signal  generator 
operating  at  576Khz.  Quadrature  (i.e.  the  two  signals  are 
out  of  phase  by  90  degrees)  signals  at  3Khz  drive  the 
elevation  and  azimuth  galvanometers. 

An  EMI7000  tape  recorder  was  used  to  record  all  test 
data.  Each  signal  output  from  the  HP467A’s  was  applied  to 
three  tape  recorder  channels  and  one  channel  of  the 
ADC/Computer  Interface.  The  gain  settings  of  the  three  tape 
recorder  channels  were  adjusted  to  accommodate  signals  that 
were  the  expected  nominal  and  plus  and  minus  50%  from 
nominal.  Other  sources  of  data  recorded  by  the  EMI7000  are 
shown  in  Figure  2. 

The  direct  memory  access,  DMA,  of  the  HP9836  computer 
was  used  as  a  transient  waveform  recorder  to  collect  real 
time  range  and  calibration  data  for  immediate  and  future  data 
reduction.  All  four  channels  of  the  ADC/Computer  Interface 
are  sampled  synchronously  at  a  36Khz  rate  and  data 
transferred  to  the  computer  by  DMA  at  144  kilobytes 


Data  Acquisition  Electronics 


Initiation  of  range  data  collection  was  synchronized  to  the 
first  positive  crossing  of  the  3Khz  signal  driving  the  red 
galvanometer  after  a  range  trigger  is  received.  This 
synchronization  was  incorporated  to  facilitate  the  phase 
measurement  required  in  processing  the  data.  When  range  data 
was  collected,  the  red,  blue,  and  3Khz  reference  signals  were 
sampled  and  stored  by  the  computer. 

Each  hologram/retroref lector  combination  was  calibrated 
by  placing  it  in  a  test  fixture  which  allowed  it  to  be 
rotated  in  two  orthogonal  planes,  simulating  pitching  and 
yawing  motion.  During  this  procedure  the  red  and  blue 
signals  and  the  x  and  y  positions  of  the  gimbal  positioning 
device  were  also  sampled  and  stored  by  the  computer. 
Appendix  A  contains  the  schematics  for  the  data  acquisition 
electronics . 


DATA  PROCESSING 


Since  data  processing  is  dependent  on  the 
characteristics  of  the  holographic  grating  and  corner 
reflector  combination,  the  reflector  assembly  will  be 
described  first. 


The  reflector  assembly,  Figure  3,  consists  of  a 
holographic  grating  and  a  cube  corner  reflector.  The 
reflector  assembly  is  small  and  rugged  with  typical  dimension 
of  6.4  mm  diameter  and  5  mm  high.  The  input  light  beam  is 
diffracted  by  the  hologram  into  the  zero  and  higher  order 
beams.  For  clarity,  only  the  zero  order  and  one  of  the 
signed  first  order  beams  are  shown  in  the  figure.  These  two 
beams  are  retroref lected  by  the  cube  corner  and  again 
diffracted  by  the  hologram.  Note  that  the  zero  and  first 
order  beams  are  again  diffracted  into  zero  order  and  first 
order  beams.  The  emerging  composite  zero  order  consists  of 
a  number  of  parallel  component  beams.  These  component  beams 
follow  different  optical  paths,  consequently  they  interfere 
constructively  or  destructively,  depending  on  relative  phase 
difference.  Figure  4  shows  the  modulation  of  the  reflection 
coefficient  due  to  optical  interference  as  a  function  of  the 
angle  of  incidence. 

The  resulting  fringe  angle  of  the  hologram  corner  cube 
combination  is  dependent  on  the  fringe  angle  of  the  hologram 
when  it  is  produced,  the  index  of  refraction  and  height  of 
the  corner  cube  reflector.  Figure  5  is  a  diagram  showing  the 
optical  set  up  for  producing  the  hologram.  The  diffraction 
angle  for  the  hologram  is 

Qd  =  2  sin-1 (X/2s) 

where  X  is  the  wavelength  of  the  laser  source  and  s  is  the 
grating  period. 

The  fringe  angle  for  the  hologram  and  corner  cube 
reflector  combination  is  approximately 

©f  =  nX/2h0d 

where  h  is  the  height  and  n  is  the  index  of  refraction  of  the 
cube  corner  reflector. 


If  the  change  in  angle  of  the  retroref lector  is 
unidirectional  then  the  amplitude  of  the  reflection 
coefficient  could  be  used  to  measure  the  angle  of  the 
retroref lector .  In  reality  the  amplitude  does  change  due  to 
effects  other  than  changes  in  the  angle  of  the  cube  corner 
reflector  (e.g.  noise,  changes  in  laser  power  output).  To 
determine  the  direction  of  change  of  the  retroref lector 
angle,  quadrature  signals  are  required.  With  quadrature 
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Figure  3.  Reflector  Assembly. 
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Figure  4.  Reflection  Coefficient. 
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Figure  5.  Technique  used  in  making  a  1°  Diffraction  Angle  Hologram. 


(i.e.  signals  are  out  of  phase  by  90° )  signals  phase  can  be 
measured,  direction  of  phase  change  determined,  and  the 
amplitude  effects  which  are  not  related  to  angle  can  be 
eliminated . 

The  spin  of  a  spinning  projectile  provides  the 
modulation  which  enables  data  reduction.  For  a  slow  spin  or 
non-spinning  projectile  the  laser  beam  is  dithered  (angle 
modulated)  in  order  to  change  the  incident  angle  of  the  input 
light  beam  at  the  retroref lector .  This  is  accomplished  by 
the  azimuth  and  elevation  galvanometers  and  beam  directing 
optics . 

Angle  modulation  of  the  beam  across  the  fringe  pattern 
shown  in  Figure  4  results  in  quadrature  relationships 
existing  between  frequency  components  of  the  photodetector 
current  which  are  harmonicly  related  to  the  modulating 
frequency.  Figures  6  through  8  support  the  analytical 
description  to  follow  from  which  the  software  algorithm  was 
developed.  This  analytical  description  is  for  one  axis  only 
and  is  not  as  rigorous  as  the  description  for  the  analytic 
model  which  was  developed  at  the  end  of  the  program. 

Figure  6a  represents  the  photodetector  output 
(equivalent  to  the  reflection  coefficient)  and  Figure  6b  the 
angle  modulation  at  the  dither  frequency.  Figure  6c  is  the 
photocurrent  for  the  angle  modulation  shown  in  Figure  6b.  I 
is  the  photodetector  current  output  as  a  function  of  the 
angle  of  the  corner  cube  reflector  and  9r  =  211  (yaw 

angle)/(fringe  angle).  Ip  is  one  half  the  peak-to-peak  value 
of  the  photodetector  current  and  t  is  time.  mp  is  the  peak 
deviation  of  the  angle  modulation  in  radians  and  wm  the 
modulation  frequency  of  the  dither  in  radians  per  second. 

The  equations  that  describe  the  modulation  effects  of 
the  dither  are 

I  ( 0r  ,  t )  =  Ip  sin(9r  +  mp  sin  wmt) 

=  Ip  sin  ©r  cos(mp  sin  wmt)  + 

Ip  cos  9r  sin(mp  sin  w*t) 

=  Ip  sin  ©r  [Jo(mp)  +  2J2(mp)cos  2  wmt  + 

2J4 (mp )cos  4wm t  +  •  .  . ]  + 

Ipcos  ©r  [2Ji(mp)sin  w»t  +  2J»(mp)sin  3w»t  +  .  .  .] 

where  the  values  of  Jn(mp)  are  determined  from  the  graph  of 
Bessel  Functions,  Figure  7.  Note  that  the  amplitudes  of  the 
odd  and  even  harmonics  above  are  quadrature  functions  The 
amplitude  of  the  odd  harmonics  is  proportioned  to  the  cosine 
of  0r  and  the  amplitude  of  the  even  harmonics  is  proportional 


b!  Phase  Deviation  (radians) 


igure  6.  Angle  Modulation  of  the  Fringe  Angle. 


Figure  7.  Bessel  Functions  for  the  first  8  orders 


to  the  sine  of  the  0r  .  The  fundamental  and  second  harmonic 
were  used  in  the  breadboard  system. 

Figure  8a  is  the  photodetector  output  with  a  graphic 
representation  showing  the  relation  between  the  fundamental 
and  second  harmonic  components  of  a  reflected  beam  for  a 
small  value  of  mp .  The  figure  shows  the  relationship  between 
the  reflected  signal  and  the  dither  signal,  when  9r  =  0,  K, 
2TT,  ■  •  ■  .  The  fundamental  frequency  dominates  the  signal  when 
6r  =0,  IT,  21T,--*  .  The  second  harmonic  dominates  the  signal 

when  0r  =  IT/2,  3TT/2,  51T/2  •••  .  Note  the  Pi  radian  phase 

reversal  between  the  positive  and  negative  slopes  for  the 
fundamental  and  the  positive  and  negative  peaks  for  the 
second  harmonic.  Figure  8b  and  c  are  the  fundamental  and 
second  harmonic  filtered  from  the  photodetector  current  by 
bandpass  filters.  The  amplitude  factor  mp  is  chosen  equal  to 
approximately  .27  resulting  in  equal  peak-to-peak  amplitudes 
for  the  fundamental  and  second  harmonic.  It  is  important  to 
note  that  the  waveforms  of  the  fundamental  and  second 
harmonic  are  in  quadrature. 

Figure  8d  shows  the  desired  quadrature  terms  as  a 
function  of  cube  corner  reflector  angle.  To  arrive  at  these 

outputs  a  phase  comparison  must  be  made  between  the 

fundamental  and  the  modulating  irequency  for  the  fundamental, 
and  between  the  second  harmonic  and  twice  the  modulating 
frequency  for  the  second  harmonic.  Devices  such  as  balanced 
mixers  used  as  a  phase  detector,  followed  by  a  suitable  low 
pass  filter,  would  produce  the  desired  output.  For  the 
computer  algorithm  the  envelope  amplitude  and  sign  is 
determined  for  each  cycle  of  reference  phase. 

After  amplitude  and  phase  detection  the  quadrature 
signals  of  Figure  8c  and  d  can  be  written  as 

El  (©r  )  =  Ip  (  2Jl  (mp  )  )C05  ©r 

E2(6r)  =  Ip(2J2(mp))sin  9r 

where  Ei  (©r  )  is  the  detected  envelope  of  the  fundamental  and 
E2  (6r  )  the  detected  envelope  of  the  2nd  harmonic.  Then 

E2  (Or  )/Ei  (6r  )  =  Ip  2J2  ( mp  )sin  ©r  /Ip  2Ji  (mp  )cos  ©r 


=  J2  sin  ©r  / J l  cos  ©r  =  (J2/J1)  tan  9r 

and 

©r  =  tan'1  [  Jl  E2  (©r  )/J2El  (©r  )  ] 

The  sign  of  the  argument  of  the  arctan  identifies  the 
quadrant  for  each  quarter- f ringe  cycle.  Comparison  of  the 
signs  of  the  numerator  and  denominator  of  the  present  and 
previous  measurements  are  used  to  identify  quadrant  boundary 
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Figure  8 


Modulated  and  Detected  Waveforms 
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crossings.  Consequently,  fractional  fringe  quarter  cycles 
can  be  measured,  and  the  boundary  direction  crossing  used  to 
sum  quarter  fringe  cycles. 
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SOFTWARE  DEVELOPMENT 


Software  was  developed  using  the  preceding  analytical 
description.  It  is  important  to  note  that  the  software  was  a 
development  effort  intended  for  diagnostics  and  for  the 
future  determination  of  the  requirements  for  complete  data 
reduction  by  a  HP9836  computer  or  the  range  computer. 

The  two  programs  (written  in  the  BASIC  language)  are 
included  as  Appendix  B.  One  software  program  is  used  for 
collecting  and  processing  range  data  when  a  range  trigger  is 
received  and  the  other  is  used  to  collect  calibration  data 
for  a  given  hologram  cube  corner  reflector  combination.  The 
calibration  program  was  derived  from  the  range  program, 
consequently  parts  of  the  program  are  not  integral  to  the 
calibration  process.  Figure  9  is  the  flow  chart  for  the 
range  data  collection  and  processing  program. 

Range  data  consisted  of  the  outputs  from  the  two 
photomultipliers  and  a  3Khz  reference  (the  3Khz  signal 
driving  the  elevation  galvanometer).  For  each  cycle  of  the 
3Khz  reference  there  are  twelve  data  samples  for  each  channel 
of  the  ADC/Computer  Interface.  Transfer  of  data  was  done 
using  the  Basic  TRANSFER  command.  Data  was  transferred  to  a 
DMA  Buffer  in  real  time  for  a  predetermined  time  after 
receiving  a  range  trigger.  At  the  end  of  this  predetermined 
time,  determined  by  the  DMA  Buffer  size,  the  binary  data  was 
converted  to  decimal  and  placed  in  a  program  array  for  use  by 
the  program. 
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Figure  9.  Computer  Program 
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Figure  9.  Computer  Program  Flow  Chart  (continued) 
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Figure  9.  Computer  Program  Flow  Chart  (continued) 
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TEST  RESULTS 


The  test  configuration  and  results  of  tests  performed  jointly  by  BAC  and 
BRL  personnel  are  summarized  in  Table  1  below.  Following  the  table  is  a 
discussion  of  the  reduced  data  listed  in  the  order  of  those  rounds  which 
produced  reducible  data. 


TABLE  1  -  ROUND  FIRING  SUMMARY 


1!  N0" 

ROUND 

TYPE 

DATA 

SET 

CONFIGURATION 

RESULTS 

K  17450 

30  mm 
spinning 

Scotchlite  on  nose, 
to  check  for 
obturation 

No  range  trigger 

ip  17452 

30  mm 
spinning 

Run-1 

Single  1°  hologram 

Red  beam  set  for  yaw 

Blue  beam  set  for  roll 
Flight  path  30  m 

Data  on  tape 

No  range  trigger 

P  17453 

30  mm 
spinning 

Run-2 

II 

No  data 

No  range  trigger 

H 

Run-3 

No  retro,  range 
trigger  check 

P 

Run-4 

No  retro,  range 
trigger  check 

n 

Run-5 

No  retro,  range 
trigger  check 

KS  17458 

30  mm 
spinning 

Run-6 

Same  as  round  17452 

Early  range 
trigger 

No  data 

m  17460 

Run-7 

No  retro,  range 
trigger  check 

20  mm 
no  spin 

Run-8 

2-axis,  1°  hologram 

Cal- 3A  Flight  path 

15  m 

Data  on  tape 
and  in  computer 

1,461 

20  mm 
no  spin 

Run-9 

Cal  -4 

II 

II 

Run-10 

No  retro,  range 
trigger  check 

11 

Run-11 

No  retro,  range 
trigger  check 
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TABLE  1  -  ROUND  FIRING  SUMMARY 


NO. 

ROUND 

TYPE 

DATA 

SET 

COff  IGURATION 

RESULTS 

17464 

30  mm 

Run-12 
Cal -5 

Set  up  for  in-bore 
measurement.  Folding 
mirror  in  blast  room 

No  data  due  to 
obturation 

17465 

25  mm 

Run-13 
Cal -6 

Muzzle  angle  and 
linear  displacement 

Collar  moved 

No  data 

setup.  Reflector  in 
collar  on  muzzle. 


17466 

25  mm 

Run-14 

Cal -6 

(1 

Displacement  and 
angle  data 

Gun  moved  back 
16mm  in  last 
two  shots 

17467 

25  mm 

Run-15 

Cal -6 

II 

Muzzle  pressure 
gauge  failed 

17468 

25  mm 

Run-16 

II 

17469 

25  mm 

Run-17 

II 

23 


Round  17452  -  The  purpose  of  this  round  was  to  test  the 
feasibility  of  determining  pitch,  roll,  and  yaw  angles  of  a 
spinning  projectile  with  a  fixed  beam  (no  dither).  The  30  mm 
round  was  fitted  with  a  single  10  hologram  and  cube  corner 
reflector.  The  calculated  fringe  angle  is  6  milliradians , 
based  on  the  hologram  diffraction  angle  and  the  specified 
height  of  the  retroref lector . 

A  retroref lector  whose  far-field  pattern  consisted  of 
two  spots  was  used  to  facilitate  roll  measurement.  The 
optics  for  the  blue  beam  were  adjusted  to  form  the  far-field 
pattern  of  the  return  beam  on  a  slit  in  front  of  the 
photomultiplier  so  an  electrical  pulse  occurs  each  time  the 
projectiles  rolls  180° .  The  reflection  assembly  on  a 
spinning  projectile  generates  a  frequency  modulated  (fm) 
signal,  where  the  number  of  cycles  per  half  revolution  is 
proportional  to  the  total  yaw  angle  (Reference  2).  The  polar 
angle  is  proportional  to  the  time  between  the  minimum 
frequency  point  of  the  fm  signal  (red  beam)  and  the 
electrical  pulse  from  the  blue  beam. 

Figure  10  shows  the  laser  measurement  results  of  total 
yaw  angle  as  a  function  of  time.  Data  obtained  from  the 
analysis  of  a  series  of  orthogonal  spark  shadowgraphs 
obtained  for  this  round  are  shown  in  Figure  11.  The  laser 
measurement  is  quite  similar  in  shape  to  the  smooth  curve 
fitted  to  the  shadowgraph  data.  It  is  interesting  that  the 
minimum  yaw  points  at  10m  are  not  on  the  smooth  curve,  but 
they  show  up  on  both  the  raw  shadowgraph  data  and  the  reduced 
laser  data. 

Figures  12a  and  b  are  laser  and  shadowgraph  measurements 
of  the  horizontal  and  vertical  yaw  angle  components.  The 
square  boxes  on  the  laser  curve  mark  the  points  at  which  the 
pair  of  far-field  blue  spots  line  up  with  the  slit.  The 
projectile  rolls  180°  between  each  pair  of  square  boxes,  thus 
the  laser  angle  sensor  can  measure  all  three  projectile 
angles . 


Round  17453  and  17458  -  No  data  was  acquired  on  these 
rounds  due  to  problems  with  the  range  trigger.  However,  the 
rounds  were  calibrated  prior  to  each  firing.  The  calibration 
technique  is  described  below. 

The  round  is  mounted  in  a  precision  two-axis  motorized 
optical  mount.  The  mount  has  a  digital  encoder  that  is 

coupled  to  the  HP-9836  computer.  The  initial  angle  is 

adjusted  so  that  the  Fresnel  reflection  from  the  reflector 

assembly  returns  to  the  laser.  The  projectile  is  rolled  so 

that  motion  on  one  axis  produces  no  modulation  of  the  fringes 
for  the  orthogonal  axis.  The  mount  is  then  driven  slowly  and 
the  fringe  pattern  versus  encoder  counts  is  stored  in  the 
computer.  Figured  13  and  14  show  typical  fringe  patterns. 
The  number  of  counts  per  fringe  is  found  by  the  geometric 
construction  shown  on  the  figure. 


.  -  571  COUNTS  - 

Red  Data  Cal  1  (counts/100) 

1240  COUNTS  -  10'2  rad 

FRINGE  ANGLE  -  4.61  mrad  (0.264  deg.) 

Figure  13.  Fringe  Angle  Calibration,  Red. 


Rounds  17460  and  17462  -  The  purpose  of  these  rounds  was 
to  evaluate  the  use  of  circular  dither  with  two  laser  lines 
to  measure  all  three  angles  of  a  non-spinning  projectile. 
The  reflector  assemblies  were  calibrated  in  the  precision 
2-axis  mount  prior  to  firing  the  rounds.  The  hologram 
configuration  was  the  orthogonal  nested  holographic  grating 
as  described  in  the  preceding  section.  The  hologram 
diffraction  angle  was  one  degree. 

Fringe  data  was  obtained  for  both  of  these  rounds,  and 
the  signals  are  quite  clean.  Unfortunately,  the  rounds  were 
very  dynamic.  The  fringe  rate  is  too  high  for  processing 
with  the  algorithm  that  was  developed  on  this  program.  The 
date  reduction  algorithm  is  designed  for  the  case  where  the 
fringe  rate  is  small  compared  to  the  dither  rate.  Figure 
15  shows  the  output  of  the  data  reduction  algorithm  for  these 
rounds.  Figure  16  is  the  range  data  for  round  17460. 


Computer  Output,  Round 


Round  17464  -  In-Bore.  The  optical  unit  was  in  the 
loading  room.  The  laser  beam  passed  through  a  50  ram  hole  in 
the  blast  room  wall  and  was  directed  by  a  pair  of  mirrors 
into  the  muzzle.  The  projectile  base  was  fitted  with  a 
rubber  seal  in  an  attempt  to  prevent  the  passage  of 
combustion  gases.  The  seal  failed,  and  no  data  was  obtained 
from  this  round  due  to  obscuration  of  the  optical  beam  by 
leakage  of  hot  combustion  gases. 

Rounds  17466  and  17467  -  Muzzle  Angle  and  Displacement. 
This  setup  is  similar  to  the  in-bore  setup,  except  that  the 
second  folding  mirror  directed  the  laser  beams  to  a  reflector 
assembly  mounted  on  a  collar  attached  to  the  muzzle.  In 
addition,  the  Boeing  Dynamic  Displacement  Measurement  System 
(DDMS)  was  used  to  measure  the  horizontal  displacement  of  the 
muzzle . 

Figures  17  and  18  are  reduced  angle  data  for  these 
rounds.  The  blue  laser  measures  the  azimuth  angle  and  the 
red  laser  measures  the  elevation  angle. 

DDMS  counter  output  versus  sample  number  is  shown  in 
Figure  19.  Positive  numbers  correspond  to  displacement 
toward  the  right  side  of  the  range,  looking  down-range.  The 
sample  rate  is  36  kHz. 

The  data  for  run  14,  round  17466,  indicates  no  motion 
for  the  first  5.6  milliseconds  (200  samples),  then  a  linear 
displacement  of  .280  millimeters  in  the  following  2.5 
milliseconds.  The  DDMS  record  then  becomes  flat,  indicating 
that  hot  muzzle  gases  have  passed  through  the  beam.  The  hot 
gas  destroys  the  lateral  coherence  of  the  signal  beam,  so  the 
fringe  count  stops.  The  hot  gas  also  affects  the  angle 
measurement  system. 

The  displacement  data  can  be  used  to  estimate  angular 
motion.  Assume  that  the  gun-tube  bending  follows  a  circular 
arc,  then  the  angular  change  at  the  end  of  the  tube  is  equal 
to  the  displacement  divided  by  one-half  the  tube  length. 
This  results  in  an  angle  of  0.56  milliradian  (0.032  degrees) 
for  round  17466. 

Proximity  probes  (used  by  BRL  personnel  following  the 
above  series  of  tests)  registered  a  peak  pointing  angle  of 
roughly  0.1  degree  in  the  horizontal  direction  and  a  somewhat 
larger  angle  in  the  vertical  direction. 


msec 


Run  15 


Muzzle  Angle,  Round  1746 


ANALYTIC  MODEL 


The  equations  relating  the  projectile  angular  motion 
with  the  output  current  from  the  photodetector  will  be 
developed  in  this  section.  The  relationships  for  the  general 
case  (a  spinning,  yawing  projectile  observed  by  an  optical 
system  employing  circular  dither)  are  first  derived.  The 
simplifications  that  result  when  using  a  fixed  laser  beam  or 
non-spinning  projectile  are  then  discussed. 

A  right-handed  coordinate  system  is  used  in  the 
analysis,  with  the  X  axis  horizontal,  the  Y  axis  vertical  and 
the  Z  axis  along  the  range  centerline.  The  optical 
centerline  also  coincides  with  the  Z  axis.  The  projectile 
undergoes  pitching  and  yawing  motions  as  a  function  of  time, 
so  the  orientation  of  a  projectile  based  coordinate  system 
(x,y,z)  may  be  specified  by  the  direction  cosines  (cosap, 
cosBp )  as  shown  in  Figure  20~a. 

The  intensity  of  light  reflected  by  the  holographic 
grating/retroref lector  combination  placed  at  the  nose  of  the 
projectile  is  a  function  of  the  angular  orientation  of  the 
input  laser  beam  relative  to  an  axis  normal  to  the  hologram 
(Reference  1) .  It  is  therefore  convenient  to  express  the 
orientation  of  the  body  in  terms  of  the  angular  coordinates 

( Up , Vp )  =  (cosOp , cosfip ) 

where  up  and  vp  may  be  thought  of  as  the  pitch  and  yaw  of  the 
projectile.  The  angles  are  chosen  so  that  up  is  oriented 
perpendicular  to  the  hologram  grating  lines  (i.e.  parallel  to 
the  direction  in  which  light  is  diffracted).  If  the  input 
angle  is  small,  the  reflection  coefficient,  R[up(t)]  is 
approximately  a  periodic  function  of  the  input  angle  as  shown 
in  Figure  20~b. 

As  part  of  the  optical  system,  circular  dither  will  be 
imparted  to  the  laser  light.  Thus  at  any  instant  of  time  the 
angle  between  a  ray  representing  the  portion  of  the  laser 
beam  striking  the  retroref lector  and  the  coordinate  axes  can 
be  defined  by  the  angular  coordinates 

(ui  ,vi  )  =  (cosoi  ,cosou  ) 

where  a  1  and  0  i  are  the  angles  between  the  incident  light  and 
the  X  and  Y  axes  respectively .  Employing  a  trigonometric 
identity,  the  angle  between  the  laser  beam  (or  the  body  axis) 
and  the  Z  axis  may  be  expressed  as 

y  =  COS'Ml  -  COS2a  -  COS20)l/2 

Figure  20-c  illustrates  the  optical  and  projectile 
coordinate  systems.  The  axes  ui  and  vi  are  centered  on  the 


optical  centerline  with  ui  horizontal.  The  body  fixed 
projectile  axes,  up  and  vp ,  rotate  as  the  projectile  rolls. 
The  circle  represents  a  circular  dither  of  the  laser  beam 
The  polar  coordinates  of  the  center  of  the  laser  beam  are 
(  n  ,6i  )  . 

The  circular  dither  can  be  resolved  into  time  varying 
components  in  the  projectile  coordinate  system 

Up  ( t )  =  rd  COS  (Wdt  -  wrt  ♦  04  o  -ir  o  )  +  Uo  ( t ) 

vp(t)  =  rd  sin  (wdt  -  wi  t  ♦  0do  -iro)  +  uo  ( t ) 

where  rd  =  angular  radius  of  the  circular  dither 

wd  =  dither  angular  frequency 
wr  =  roll  angular  frequency 
0do  =  initial  dither  angle 
0r o  -  initial  roll  angle 

uo  ( t )  ,  \y(  t )  =  position  of  optical  axis  in  projectile 
coordinates . 

The  component  of  the  dither  in  the  up  direction 
(perpendicular  to  the  holographic  grating  lines)  produces  a 
modulated  output,  while  the  component  in  the  vp  direction 
(parallel  to  the  grating  lines)  does  not 

The  detector  current  is  the  product  of  two  time  varying 
functions.  Since  the  intensity  of  the  laser  is  not  spatially 
uniform,  as  the  beam  is  swept  across  the  reflector  (due 
either  to  the  dither  or  the  motion  of  the  projectile)  the 
light  entering  the  retroref lector  is  amplitude  modulated. 
The  reflector  further  modulates  the  beam  because  the 
reflection  coefficient,  R[up(t)],  varies  with  the  angle  of 
incidence . 

If  it  is  initially  assumed  that  the  laser  beam  profile 
is  not  a  function  of  z,  then  th<-  intensity  may  be  represented 
by 

I  ( t )  =  (  rb  ,  0b  ) 

where  rb  =  ( rd 2  +  n  2  -  2rdricos(0do  ♦  wdt  -  01  ))i/2 

0b  =  cos-i[(rd2  -  rb2  -  n2)/2rbri] 

The  circular  dither  is  actually  a  conical  scan  with  a 
fixed  point  at  a  distance  z  from  the  initial  projectile 
position.  The  linear  displacement  of  the  beam  with  respect  to 
the  retroref lector  is  a  function  of  z-zo  (when  z-zo  there  is 
no  translation  of  the  beam).  The  z  axis  dependence  can  be 
placed  into  the  expression  for  I ( t )  by  replacing  n  with 
n  ( t ) ,  where 

n  ( t )  -  ri  (  0  )  •  (  zo  '  z  ( t )  )  /zo 
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The  photodector  current  is  therefore 


i  ( t )  =  Pri(t)  R[ up ( t )  ] 

where  P-  detector  responsivity  (amp/watt) 

T -  optical  efficiency 
I  =  laser  beam  intensity 

Two  special  cases  will  be  considered,  both  assuming  that 
the  laser  beam  is  uniform,  so  that  I(t)  is  constant. 

(1)  The  retroref lector  reflection  coefficient  is  a  pure 
sinusoid  with  period  8f  . 


(  2  ) 


The  reflection  coefficient 
higher  order  components. 


is  periodic,  but  has 


Case  1:  R  [  up  ( t  )  ]  -  a  ♦  b  s in {  21fup  ( t ) /8f  -  c) 


where  a  =  average  value 

b  =  amplitude 
0f  =  angle  between  fringes 

c  =  phase  at  t-0 

up  i  t  )  -  rd  cosfwdt  -  wr  t  +  9do  -  8r  o  )  +  uo  ( t  ) 

Let  d  -  Pn,  k  =  2  IT  /0f  and  w  =  wd  t  -  wr  t  +  0d  o  -  0r  o 

then  the  detector  current  is 


i(t)  =  ad  ♦  bd  sin[krdcos  w  +  kuo  ( t )  -  c] 

Express  the  sine  term  as  a  sum  to  get 

i  t  -  ad  +  bd  sin ( krd  cos  w )  cos ( kuo ( t )  -  c ) 

+  bd  cos(krdcos  wj  sin(kuo(t)  -  c) 

We  can  then  expand  the  terms  with  sinusoidal  argument  using 
the  Bessel  function  identities 


sinf  zcosei  -  2Ji(z)cos0  -  2J3(z)cos30  +  ... 
cos ( zcos0 1  =  Jo ( z  )  -  2J2(z)cos20  +  2J<(z)cos40  - 
The  first  few  terms  are 


i(t)  =  ad  ♦  bd  Jo ( krd  )  sin(kuo(t)  -  c) 

♦  2bd  Ji(krd)  cos  w  cos(  kuo  ( t )  -  c  ) 

-  2bd  J2(krd)  cos  2w  sin(kuolt)  -  c) 


2bd  Jsikrd)  cos  3w  cos(kuo(t 


c  *  ♦ 


These  harmonic  components  can  be  separated  by  filtering,  then 
bidirectional  changes  in  uo(t)  can  be  determined  by  tracking 
the  coefficients  of  cos  [kuo(t)]  and  sin  [kuo(t)]  for  a  pair 
of  even  and  odd  harmonics .  Note  that  the  harmonic 
frequencies  are  multiples  of  the  difference  between  the 
dither  and  roll  frequencies.  If  the  dither  is  high  compared 
to  the  roll  rate,  roll  produces  a  slow  change  in  the  phase  of 
the  dither. 

Case  2:  R(u)  is  periodic  but  not  a  sinusoid. 

The  projectile  reflector  must  be  calibrated  before 
firing.  During  calibration,  R(u)  is  measured  over  at  least 
one  period  of  6f  .  If  the  dither  is  on,  the  amplitudes  of  two 
adjacent  harmonics  are  stored  in  a  look-up  table  of  harmonic 
amplitudes  versus  angle.  If  the  dither  is  off,  the  table 
entries  can  be  calculated  as  follows: 

Write  R(u)  as  a  Fourier  series 

R(u)  =  ao /2  +  ai  cos(ku)  +  a2  cos(2ku)  +  ... 

+  bi  sin(ku)  +  b2  sin(2ku)  +  ... 

Replace  u  with  u  +  rd  cos  (w),  then  expand  each  term  using  the 
Bessel  function  identities  and  collect  like  harmonics.  The 
first  few  term  of  the  expansion  and  rearrangement  are: 

m 

Ro(u)  =  ao/2  +  E  [aj  cos(jku)  +  bj  sin(jku)]  Jo  ( jkrd  ) 

j  =  l 

m 

Ri  (u)  =  2  cos  w  £  [aj  sin(jku)  -  bj  cos(jku)]  Ji(jkrd) 

j  =  l 

m 

R2(u)  =  2  cos  2w  E  [aj  cos(jku)  +  bj  sin(jku)]  J2(jkrd) 

j  =  l 

m 

R3(u)  =  2  cos  3w  £  [aj  sin(jku)  -  bj  sin(jku)]  J3(jkrd) 

j  =  l 

These  are  the  harmonic  amplitudes  for  the  look-up  table.  The 
table  is  used  to  reduce  the  data  after  the  round  is  fired. 

The  data  is  reduced  in  increments  of  the  dither  minus  roll 
period.  For  each  increment,  harmonic  component  amplitudes 
can  be  found  with  the  digital  filtering  routine  in  the  data 
reduction  algorithm.  The  look-up  table  is  used  to  find  the 
fractional  angle  u,  where  u  £  9f  .  The  total  angle  is 
incremented  or  decremented  by  0 t  each  time  u  moves  through 
the  table. 


PRELIMINARY  DESIGN 


This  section  contains  preliminary  design  information  for 
a  three  angle  measurement  system.  The  design  addresses  the 
limitations  of  the  breadboard  system  tested  at  BRL,  and 
defines  the  areas  which  require  further  development  before  a 
final  design  can  be  implemented. 

The  major  limitation  of  the  breadboard  system  was  the 
dither  rate.  An  available  galvanometer  deflection  system  was 
used  to  generate  a  conical  scan  at  3  KHz.  However,  in  many 
cases  the  fringe  rate  was  comparable  to  the  dither  rate.  The 
dither  rate  must  be  much  higher  than  the  fringe  rate, 
otherwise,  the  desired  effect  of  the  dither  modulation  is 
destroyed  and  angle  information  cannot  be  recovered.  Spin 
stabilized  rounds  with  fringe  rates  of  30  KHz  have  been 
observed  in  our  test  program.  A  dither  rate  of  the  order  of 
1  MHz  would  be  appropriate.  Oscillating  mirrors 

(galvanometer  type  of  resonant  scanners)  have  a  capability  up 
to  20  KHz.  Polygonal  mirrors  on  turbine  motors  can  generate 
higher  scan  rates,  but  produce  linear  rather  than  conical 
scans.  The  mechanical  limit  (reference  3)  for  a  polygonal 
mirror  is: 

w  *  - f  8UTS 

21Tr°  yp(3+n) 

where  ro  =  Distance  from  center  to  edge 
UTS  =  Ultimate  tensile  strength 
n  =  Poisson’s  ratio  of  material 
p  =  Volumetric  density 

Substituting  the  values  for  6061  T6  aluminum,  the  mechanical 
limit  for  a  5  cm  diameter  scanner  is  about  104 
revolutions/second.  If  the  mirror  had  20  facets,  the  scan 
rate  would  be  200  KHz.  The  practical  limit,  due  to  mirror 
deformation  and  safety,  is  well  below  the  mechanical  limit. 

Acousto-optic  beam  deflectors  (AOD)  are  an  alternative 
method  for  generating  a  conical  scan.  Two-axis  devices  are 
available,  so  a  circular  dither  is  possible.  The  limit  is 
set  by  the  transit  time  of  a  sound  wave  across  the  aperture 
of  the  device.  The  transit  time  is  the  amount  of  time 
required  to  move  the  deflected  beam  from  one  position  to 
another.  Typical  valves  for  commercial  deflectors  are  of  the 
order  of  10  microseconds .  In  the  breadboard  system  the 
return  beam  was  sampled  at  12  positions  around  the  dither 
circle.  If  an  AOD  is  used  with  10  microsecond  transit  time 
to  sequentially  deflect  the  beam  to  12  points  on  a  circular 
arc,  the  maximum  rate  would  be  8.51  KHz.  This  rate  is  also 
inadequate  for  a  spinning  projectile. 


A  solution  to  the  high  dither  rate  requirement  is  to  use 
a  circular  array  of  diode  laser  beams,  where  the  diode  lasers 
are  sequentially  modulated  to  produce  a  rotating  optical 
source.  Very  high  dither  frequencies  are  possible  with  this 
technique.  (We  are  currently  modulating  a  Mitsubishi  model 
ML  3001  diode  laser  at  a  frequency  of  1  GHz.)  The  wavelength 
of  laser-diodes  are  temperature  sensitive.  Diode  lasers  are 
available  with  fiber  optic  pigtails,  so  the  lasers  can  be 
mounted  on  a  temperature  controlled  heat  sink  and  the  fibers 
led  out  to  form  a  circular  array  in  the  focal  plane  of  the 
collimating  lens.  Two  color  operation  is  accomplished  by 
selecting  lasers  so  the  even  and  odd  fibers  radiate  in 
different  wavelength  bands,  such  as  830±5nm  and  850  ±  5nm. 

Figure  21  shows  two  optical  configurations.  Both 
configurations  show  an  off-axis  paraboloid  mirror  for 
collimating  the  beam.  The  beam  diameter  would  be  of  the 
order  of  30  cm,  so  the  refracting  optics  as  used  in  this 
program  become  more  expensive  than  off-axis  optics. 
Furthermore,  the  Fresnel  reflection  from  the  center  of  a 
refractor  is  avoided  with  reflecting  optics. 

Figure  21a  is  a  three-axis  system  for  measurements  of 
spinning  projectiles,  using  techniques  verified  on  this 
program.  The  concept  shown  in  Figure  21b  is  a  circular 
dither  system  using  diode  laser,  as  described  above.  The 
first  technique  is  lower  risk  and  cost,  since  all  the 
techniques  have  been  demonstrated  in  this  program. 


PHASE  II  PROGRAM 


The  major  task  in  Phase  II  is  the  design  and  fabrication 
of  a  three-component  angular  measurement  system.  The  results 
of  the  Phase  I  test  program  made  a  significant  impact  on  the 
requirements  for  the  Phase  II  system.  The  major  requirement 
is  for  a  much  higher  dither  rate  in  order  to  follow  the 
dynamics  of  some  of  the  projectiles.  Other  requirements  are 
increased  laser  beam  diameter,  improvements  in  the  data 
reduction  algorithms,  and  increased  stability  of  the 
holograms  in  the  humid  environment  at  BRL.  The  Phase  II 
tasks  are  listed  below: 

1)  Design  and  fabricate  a  three -component  angular 

measurement  system.  Two  variations  of  the 

three-component  system  are  possible. 

a.  Total  yaw  system  for  spinning  projectiles.  This 

technique  was  demonstrated  on  the  prior  contract, 
and  improved  to  measure  all  three  angles  in  the 
present  contract. 

b.  Three-component  system  with  circular  dither  for  both 

spinning  and  non-spinning  projectiles.  The 

technical  difficulties  are  greater  and  the  risks  are 
higher  for  this  system.  At  this  time  it  appears 
that  the  only  method  for  generating  the  required 
high  dither  rate  is  to  use  a  circular  array  of 
optical  fibers,  each  driven  by  a  diode  laser,  where 
the  relative  phases  of  the  diode  lasers  is  adjusted 
to  produce  a  conical  scan.  The  dither  generation 
technique  should  be  thoroughly  evaluated  with  a 
breadboard  device  prior  to  fabrication  of  the 
three-component  system. 

2)  Data  reduction  algorithms 

a.  No  computer  data  reduction  algorithm  has  been 

developed  for  la.  Such  a  system  would  require  the 
development  of  a  suitable  algorithm. 

b.  The  algorithm  for  lb  has  partially  been  developed  in 
Phase  I .  Additional  development  is  required  based 
on  the  results  of  Phase  I  and  the  analytic  model 
developed  in  Phase  I.  One  of  the  important  problems 
to  be  addressed  is  the  magnitude  of  error  caused  by 
the  deviation  of  the  reflector  reflection 
coefficient  from  a  sinusoid.  The  algorithm  already 
developed  is  for  the  sinusoidal  case,  but  the 
reflection  coefficient  is  approximately  sinusoidal 
only  when  the  diffraction  efficiency  is  low.  In 
order  to  maximize  the  return  power,  it  is  necessary 
to  operate  with  a  reflection  coefficient  versus 
angle  that  departs  significantly  from  a  sinusoid. 
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development  would  be  a  joint  Boeing  and 


-echniques  to  stabilize  the  holograms  after 
The  diffraction  efficiency  of  the 
tanged  with  time  due  to  the  change  in 
dity  at  BRL.  Some  of  the  holograms  became 
lusable . 

item  of  lb,  design  and  fabricate  signal 
timing,  and  computer  interface 
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! BRL_6_28 

i - 

{PROGRAM  FOR  OBTAINING  THE  CALIBRATION  CURVE  FOR 
! A  GIVEN  HOLOGRAM  CORNER  CUBE  COMBINATION. 

i - 

Comments l! This  program  uses  the  difference 

{method  for  determining  the  amplitude 
{for  f  2  4  f  1  . 

!The  largest  Num_samp les  is  32K. 

{Graphing  is  included  for  200  enve- 
! lopes/graph .  a  of  graphs  depends  on 
! the  number  os  Num_samples. 

{Parameters  of  consideration:  Ratios 
•of  f2/f1,  degrees/count  red(Dpcr), 
!degrees/count  blue(Dpcb) 

!The  doppler  effect  is  included  but 
!not  used(It  was  developed  when  this 
{program  was  configured.).  It  uses 
!the  parabolic  fit  to  determine  the 
{location  of  the  peak  using  the  greater 
!of  the  ampplitudes  of  fl  4  f2.  The 
idoppler  correction  is  indepentent  for 
!the  red  4  blue. 

!For  a/d  interface-  reset  for  counter 
land  normal  cal  switch,  up  for  cal. 

PRINTER  IS  1 
PRINT  CHRSM2) 

PRINT  TABXY(30, 15), "CALIBRATION  PROGRAM- 
WAIT  3 

PRINT  TABXY (30,15),“ 

Initialize:  ! 

PRINTER  IS  1 

PRINT  CHRS(12)  '.CLEAR  SCREEN 
ON  ERROR  GOTO  Err_rout 
. OPTION  BASE  0 
GRAPHICS  OFF 
RAD 

ON  KEY  4  LABEL  -GRAPH", 15  GOSUB  Disp_out 
ENABLE 


33  Gpio'12 

84  Hpib*:704 

85  {CONSTANTS 

86  Sf *36000  {Sample  frequency 

87  Gf*3000  fgalvo  frequency  set  by  clod 

38  K1-12  {samples  at  fl 

89  K2-6  {samples  at  f2 

30  K3*4  {samples  at  f3 

91  iSianal  generator  frequency  576.000Khz 

92  {VARIABLE  LIST 

93  !Num_samp  :  number  of  samples  to  take  at  36k  z  rate 

94  ! Fb  :  filter  bandwidth 

95  !N  :  lHTCSf/fb> 

96  !Nnn  :INT  ( Num_samples/48 ) *  1 2 ,  used  in  DMA 

97  !Br(  )  :  array,  sampled  data,  red 

98  ! B 1 (  )  :  array,  sampled  data,  blue 

99  ! I f 1  :  interpolation  factor  for  fl 

100  !If2  :  interpolation  factor  for  f2 

101  !F1(  )  :  convolution  multipling  function  for  fl 

102  ! F 2 (  )  :  convolution  multipling  function  for  f2 

1 OJ  {Start  :  starting  point,  determined  by  N 

104  !Qr  :  starting  sampled  data  point  for  next  K1 
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1 05  .'interval,  including  doppler  correction,  red. 

106  !Qb  :  same  as  Or  except,  blue 

107  !Xx  :  Xx*1,  loop  for  red:Xx*2  loop  for  blue 

108  !B1(  )  :  Gimbal  position  for  sample,  red 

109  IB2(  )  :  Gimbal  position  for  sample,  blue 

110  IB1K1  :  running  sum  of  B1(  )  in  K1  interval 

111  IB2K1  :  running  sum  of  B2(  )  in  K1  interval 

112  ! F i 1 1 _ f  M  J )  :  convolution  elements  in  K1  interval,  fl 

113  ! F 1 1 1 _ f 2 ( J )  :  convolution  elements  in  K1  interval,  f2 

114  ! B 1 ( P )  :  average  of  Bit  )  in  K1  interval 

116  !B2(P)  :  average  of  B2(  )  in  k2  interval 

116  !P  :  K1  interval  of  computation 

117  !Ror  :  running  sum  of  doppler  correction,  red 

118  !Rob  :  running  sum  of  doppler  correction,  blue 

119  !Rotr(P)  :  running  sum  stored  as  a  function  of  P 

120  !Rotb(P)  :  running  sum  stored  as  a  function  of  P 

121  IQccr(P)  :  doppler  correction  for  respective  K1  interval,  red 

122  !0ccb(P)  :  doppler  correction  for  respective  K1  interval,  blue 

123  IQpsr  :  quadrant  past,  red 

124  IQpsb  :  quadrant  past,  blue 

125  IQptr  :  quadrant  present,  red 

126  IQptb  :  quadrant  present,  blue 

127  !Qsr  :  running  sum  of  quadrants,  red 

128  !Qsb  :  running  sum  of  quadrants,  blue 

129  !Tar  :  total  angle  red,  K1  interval 

130  !Tab  :  total  angle  blue,  K1  interval 

131  !Ta_red(P)  :  total  angle  red,  P  K1  interval 

'32  ! T a _ b 1 ( P )  :  total  angle  blue.  P  K1  interval 

133  IWhole  quadrant  update  table 

134  DIM  Qt (5,4) 

135  Qt ( 1 ,2 )  =  1 

136  Gt(1,4)«-1 

137  Gt(2,1)»-1 

138  Qt ( 2 , 3 )  *  1 

139  Qt(3,2)  — 1 

140  Qt  ( 3 , 4 )  =  1 

141  Qt ( 4 , 1  )  *  1 

142  Qt ( 4 , 3 )  =  -  1 

143  Qsr=10  IQuadrant  sum  red  starts  at  10 

144  Qsb=10  IQuadrant  sum  blue  starts  at  10 

145  ASSIGN  »Gp 10  TO  12;W0RD 

146  IPRINT  "  REGISTER  STATUS  BEFORE  DMA” 

147  IGOSUB  Check_stat 

148  Main:!  Main  program 

149  PRINT  CHRS (12) 

150  GQSUB  Num_samp  !»  of  samples  to  take 

161  IGOSUB  Test 

'52  IGOSUB  File_data  !  Copy  data  from  Buffer 

1  5 is  I  to  Disk  File. 

154  IGOSUB  File_array 

155  GOSUB  F 1 1 t _ bu>  int  llnput  filter  bandwidth  &  interpolation  factor 

156  BEEP 

:S7  INPUT  "If  calibration  run,  Type  C.  else,  Type  S,  for  stored", Aa$ 
158  IF  Aa$=" C"  OR  AaS-"c”  THEN 
’SS  GOSUB  Dma_transfer 

160  GOSUB  File_data 

161  GOTO  174 

162  END  IF 

>63  IF  AaS* "S"  OR  Aa$*"s"  THEN 

164  DISP  "INSERT  DATA  DISK  -  PRESS  CONT  TO  CONTINUE' 

165  PAUSE 


$4 


CAT 

GOSUB  File  array 
GOTO  174 
END  IF 
BEEP 
BEEP 
BEEP 

GOTO  157 

GOSUB  Mult  fun  !C 


GOSUB  Mul t_f un  (Generate  multiplier 

! f  unc  t ions . 

GOSUB  Comp_angle  (Convolution,  envelop  detection,  doppier  correction, 

(4  Arctan  for  red  and  blue 

DISP 

GOSUB  Disp_out 
GOTO  1 

(+•*■  +  +  +  +  +  +  +  +  +  +  +  +  *  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  4-  +  +  + 


182  Num_samp:  (Input  the  total  number  of  samples 


(to  be  taken.  Maximum  of  32K 

BEEP 

Num_samples‘ 30000 
GOTO  188 

INPUT  “TOTAL  NUMBER  OF  SAMPLES  TO  BE  TAKEN?" ,Num_samples 

ALLOCATE  REAL  Ta_red( INT (Num_samples/K1 ) ) , Ta_bl ( INT (Num_samples/K1 ) ) ,Rot ( 


NT ( Num_samples/K 1 ) ) 


RETURN 


1 90  Dma  transfer  :  ! 


DISP  “  COLLECTING  DATA" 

S  =  0 

INTEGER  Kk ( 0 : 47 )  BUFFER 
ASSIGN  SBuf  TO  BUFFER  KM*> 

(RESET  DATA  TRANSFER  ELECTRONICS 
CONTROL  Gpio,2:2  (SET  CTL0  4  CTL 1 


TO  1  & 
&  CTL  1 


CONTROL  Gpio,2;0  (SET  CTL0  TO  1  &  CTL 1  ' 

CONTROL  GPio,2;2  (SET  CTL0  4  CTL 1  TO  1 

(END  RESET 

CONTROL  §Buf,4;0 

CONTROL  Gpio,0;2 

CONTROL  GPio,3;0 

(SET  TRANSFER  LOW 

CONTROL  Gpio,2;3  (SET  CLT0  TO  0,  XFER  0 
(BEEP 

TRANSFER  SGpio  TO  @Buf 
FOR  X-0  TO  47  STEP  4 
Br ( S ) *Kk ( X ) - 1 28 
Bl(S)-Kk (X+1 HI  28 
B1 ( S )  'Kk ( X+2 )  + 1 
B2( S ) -Kk ( X+3 ) + 1 

!  PRINT  S,Br(S),BKS),B1(S),B2(S) 

S*S  +  1 
NEXT  X 

IF  S-Nnn  THEN 
GOTO  220 
ELSE 

GOTO  196 
END  IF 

CONTROL  Gpio,2:2  (SET  CLT0  4  CLT1  TO  1 
CONTROL  Gp i o , 1  ;  3 


CONTROL  Gpic 
DISP 
(PRINT  " 
(GOSUB  Check 
BEEP 


REGISTER  STATUS  AFTER  DMA" 
sta  t 


226  RETURN 

246  F i 1 t_bu_i nt :! Input  filter  bandwidth  in  hertz 

247  BEEP 

260  INPUT  “WHAT  IS  THE  DESIRED  FILTER  BANDWIDTH  IN  HERTZ?-, Fb 
2S1  N«INT(Sf/Fb) 

253  Nnn“INT(Num_samples/48)*12 

254  ALLOCATE  INTEGER  Br  ( 0  :Nnn ),  INTEGER  Bl<0:Nnn) 

255  ALLOCATE  INTEGER  B 1 ( 0 :Nnn ), INTEGER  B2(0:Nnn> 

256  ALLOCATE  INTEGER  Rot r ( 0 : Nnn ) , Ro tb( 0 :Nnn ) ,Qccr ( 0 :Nnn ) ,Qccb( 0 :Nnn ) 

25 7  PRINT  CHR$( 1 2 ) 

258  BEEP 

262  INPUT  “Real  samples  fl-12,  interpolation  factor?", I f 1 

263  BEEP 

264  INPUT  "Real  samples  f2-6,  interpolation  factor?*, If2 

265  BEEP 

266  ALLOCATE  REAL  Filt_f 1 ( 1 :Kl*If 1 ) ,Filt_f2< 1 :K1*If2) 

267  RETURN 

268  Test:!  Generate  data  -  needs  to  be  modified  for  this 

269  '.program 

272  BEEP 

273  INPUT  “What  is  the  Envelop  Peak  Ampl itude?” ,Epa 

274  BEEP 

275  DISP  "  GENERATING  DATA". 

276  Peak„value_data*0 

277  Ga l_cyc_f ra=200 

278  Ab=0 

279  ALLOCATE  INTEGER  Mmm ( 0 : Num_samples- 1 ) 

280  FOR  X-0  TO  Num_samP les- 1  STEP  4 

281  Ac*Epa*SIN( <PI*(+Ab) ) /(6*Gal_cyc_f rg )  +  (2 . 7 )  ’•SIN ( 2*PI*Ab*Gf /Sf ) ) 

282  FOR  Y-X  TO  X+3 

283  Mmm< Y ) « INT ( Ac > ♦ 1 28 

284  IPRINT  Y.Mmm(Y) 

285  OUTPUT  «Buffer  USING  “#,B" ;Mmm( Y ) 

286  NEXT  Y 

287  ...  Ab-Ab+1 

288  NEXT  X 

289  DISP 

290  GOSUB  Check_stat 

291  WAIT  5 

292  RETURN 

293  Mult_fun:  'Generation  of  the  multiplying  function. 

294  (Multiplying  function  will  be  ( SIN ( X) / ( X ) )*COS ( Y  ) 

295  ! PRINTER  IS  705 

296  (ALLOCATE  REAL  FI ( 0 : (N*If 1) )  .REAL  F2 1 0 : ( N* I f 2 )) 

297  DIM  FI (0:2000) ,F2(0:2000) 

298  RAD 

299  BEEP 

300  (Generate  1/2  multiplier  function  for  FI 

301  BEEP 

302  DISP  "Building  multiplier  function  for  FI" 

303  FI ( 0  )  - 1 

304  (PRINT  "FK0)  =  ":F1(0) 

305  FOR  P-1  TO  N *- 1 f  1 

306  Aa* (PI*P)/(N*If 1 ) 

307  Bb- ( Aa*2*Gf  ) /Fb 

308  (PRINT  “P-U:P 

309  (PRINT  "N*"  ;N*-If  1 

310  (PRINT  "(PI*P)/N-" :Aa 

3 1  1  (PRINT  "  ( Aa*2*-GF  )  /FB' " :  Bb 

312  FI ( P ) » ( SIN ( Aa ) /Aa ) *C0S ( Bb ) 

313  (PRINT  "FI ( " ;P; " ) *" ;F1 (P) 


314  NEXT  P 

315  DISP  "Building  multiplier  function  for  F2" 

316  IGenerate  1/2  of  multiplier  function  for  F2 

317  F2(0  )  -  1 

318  SPRINT  "F2(0)*";F2(0) 

319  FOR  P-1  TO  N*If2 

320  Aa= ( PI*P ) / ( N* I f 2 ) 

32 T  Bb-(Aa*4*Gf  )/Fb 

322  SPRINT  "P-“ ;P 

323  SPRINT  “N=":N 

324  SPRINT  "<PI*P)/N«":Aa 

325  SPRINT  ■ ( Aa*4*GF ) /FB* " ; Bb 

326  F2(P)=<SIN(Aa)/Aa)*C0S(Bb) 

327  SPRINT  *F2<  n ;P: " ;F2(P) 

328  NEXT  P 

329  RETURN 

330  ! 

331  ! . 

332  SConvolut ion  and  envelop  detection 

333  Sat  K1  intervals  for  fl  and  f2. 

334  ! . 

335  ! 

336  Comp_ang 1 e : ! F i 1  ter ,  envelop  detection. 


337 

Sdoppler  correction  and 

338 

Sarctan  for  red  &  blue 

34  0 

RAD 

34  1 

Tot  anqle_peak=0 

342 

P-0  !  In i  1 1 a  1  i  za ton  of  array  for 

the 

out- 

343 

Sput  anale  for  red  and  blue 

lasers . 

344 

SDetermine  first  data  point  to  start 

con- 

345 

Svolution  considering  the  value 

of  N 

and 

346 

Sthe  need  of  a  +  slope  for  fl  for  synchr- 

34  7 

Ionization. 

348 

Ouot lent* ( N+ 1 )  DIV  K1 

349 

Start* (Guot ient  +  1  )♦ 1 2 

351 

BEEP 

353 

DISP  M  C  R  U  N  C 

H  I 

N  G 

355 

PRINT  "Start*" :Start 

356 

Qq-0  {Initialize  graph  index.  Increment 

357 

! Oq  for  each  200  computations 

361 

Qr-Start+2-3 

362 

Qb*Gr 

363 

FOR  Xx*1  TO  2  !  1  *  red . 2*b 1 ue 

if  L 

IE  X >•  =  1  THEN 

365 

G*Gr 

366 

ELSE 

"67 

G*Gt> 

366 

END  IF 

369 

&  1  k  1-0 

370 

B2i 1-C 

37  3 

1 

37a 

Cor-  f  1 : ! 

7  c 

itonvoiution  over  Ft  interval 

tor 

376 

!  F  1 

377 

l 

376 

J* 

7  79 

Peak  value  f  1  *  0 

F0F  > *6  TO  -1 

■  ;  * 

;f  X  *  1  Then  pH  ’-R1  '  X  i  +BH  1 

if  x>-^  Then  B2ki*bL<X)*B2** 

383 

FOR  C-0  TO  If  1  -  1 

384  Sum_f 1 *0 

385  IF  Xx-1  THEN  ! Xx- T .Red:Xx-2.Blue 

386  FOR  M-X  TO  X+N 

387  Produc t_f 1  *Br (M+ 1 )*F 1 ( ( M-X ) *1 f 1 -C+ If  1 ) +Br I 2+X-M )*F 1 ( ( M-X ) «I f  1  +C ) 

388  Suin_f  1  *Sum_f  1 +Produc  t_  f  1 

389  !  PRINT  Q; TAB( 1 0 ) ;X; TAB(20 ) ;C; TAB(30) ;H;TAB(40) ; J; TAB(45) ; ( (M-X )« 

If  i  +C ) : 

390  !  PRINT  TAB(50>;((M-X>*If1-C+Ifn; 

391  !  PRINT  TAB160) ; (2+X-M-1 ) 

392  NEXT  M 

393  FI  9F 

394  FOR  M-X  TO  X+N 

395  Product_f1«BKM+1 )*F1 ( (M-X)«If 1-C+If 1 )+Bl(2«X-M)«Fl ( <M-X)«If  1+0 

396  Sum_f 1 -Sum_f 1 +Product_f 1 

397  !  PRINT  Q : TAB( 1 0 ) : X ; TAB ( 20 ) ; C; TAB ( 30 ) ;M; TAB ( 40 ) ; J: TAB ( 45 ) ; ( ( H-X ) 
+  If  1 +C ) ; 

398  !  PRINT  TAB( 50 ) ; ( (M-X )*If 1 -C+If 1) : 

39S  !  PRINT  TAB(60) ; (2+X-M-1 ) 

400  NEXT  fl 

401  END  IF 

402  F i 1 t_f 1 ( J) “Sum_f 1 

403  IPRINTER  IS  705 

404  'PRINT  "F 1 1  t_f 1 ( “ : J; "- ) * ;F  j 1 1  fl(J) 

405  IPRINTER  IS  1 

406  IF  ABS(Filt_f 1 (J) > >Peak_value_f 1  THEN 

407  Peak__va  lue_f  1 -ABS ( F  1 1 1  _  f  1  (J)  ) 

408  END  IF 

409  J=J+1 

410  NEXT  C 

411  NEXT  X 

412  IF  Xx*  1  THEN  B1  ( P ) - INT < B1  lc  1  /1 2 ) 

413  IF  Xx-2  THEN  B2(  P )  *  INT  ( B2H /1 2 ) 

419  ! . 

420  Env_f 1 :! Envelop  detection  over  FI 

421  .  ! i nterval  for  f 1 

422  ! . 

423  Dif f_plus_f 1 “0 

424  Dif f_minus_f 1 ”0 

425  FOR  E-Ifl+1  TO  5*If1+1 

426  Dif  f_f 1 *F 1 1 t_f 1 <E)-F ilt_f1 *  E ♦ I N  T  <  1  fi  *M /I  i  i 

427  IF  Di f f _f 1 >Di f f _p 1 us_f 1  THEN 

428  0 i f  f  _p 1 us_f 1  *  D i f  f  _f 1 

429  END  IF 

430  IF  Dif  f _F  1  < D j  f  f_runus_f  1  THEN 

431  D i f f_m i nos  f 1 *D i f f _f 1 

432  END  IF 

^33  ! PRINT  “e*";E.‘Diff_f  1*":DiM_" 

434  !PRINT  * D i f f _p 1 us  f 1  *  * : D i f f _p i us _ f ' 

635  IPRlNT  ”D i f f _m j nus  f'*“;Dif*  minus  •’ 

436  NEXT  E 

437  IValue  of  Envelop  at  P  interval  for  f' 

438  IF  D i  ♦  f  _p  1  us_f 1  > ABS t  D  i  f  f  _r. .  nu>:  *  '  T Hf  i, 

439  E nve 1 op_f 1 ■ D i f f _p 1 us_f 1 

440  ELSE 

441  E n ve 1 op_ f 1 ■ D i f f _m i nus  f' 

442  END  IF 

476  ! . 

4  79  Cor,_f  2  :  1 

460  •  6  nr.  VO  i  o  t  i  Or  Cvf  r  *  *  ,  r.  •  e  -  „  d  •  1  . 

4  -■  '  ! . 

483  j*  ■ 


63 


it". 

57 

38 

39 

30 

3’ 

32 

It  2  +  C  ) 

33 
94 
99 
5b 

37 

38 

39 
30 

1 1  2  *  C  i 

3’ 

32 

33 

34 

39 


Pealc_value_f  2*0 
FOR  X-Q  TO  Q+K1-1 
FOR  OO  TO  If 2- 1 
S  uit> _ f  2*0 

IF  Xx-1  THEN  !  Xx*  1 .Red : Xx-2 , Blue 
FOR  M-X  TO  X*N 

Product  f2-Br  (M*  1  >*F2<  ( M- X ) * I f 2-0 If 2 )  *Br  1 2«X-M)*F2(  ( M-X  >♦  If 2*C.  > 
Sum_f  2*5um_f 2*Produc  t_f 2 

!  PRINT  Q : TAB ( 1 0 ) ; X : TAB (20 ) ;C : TAB( 30 ) ;M; TAB(40 ) ; J; TAB( 45 ) ; < (M-X ) 

!  PRINT  TAB<50>;< (M-X)*If2-C*If2); 

!  PRINT  TAB(60):<2*X-M-1 ) 

NEXT  M 
ELSE 

FOR  M-X  TO  X*N 

Product  f  2-B1  (M*  1  ) *F2(  <M-X)*If2-C*If2)*Bl  (2*X-M)*F2(  (M-X)* If 2+0 
Sum_f 2-?un_f 2*Produc  t_f 2 

!  PRINT  0 : TAB ( 1 0 ) : X : TAB( 20 )  ;C: TAB ( 30 ) ;M; TAB( 40 ) ; J ! TAB( 45 ) I ( l M-X ) 

!  PRINT  TAB ( 50  >  :  (  ( M- X  >« If 2-OIf  2 L: 

'  PRINT  TAB160) ; (2*X-M-1 ) 

NEXT  M 
END  IF 

F  j j  t_f 2<  J  )  -Sum_f2 
! PRINTER  IS  70S 

! PRINT  "F  i  l t_f  2  <  *;J:"*)*:Fiit_f2(J) 

•PRINTER  IS  1 

IF  ABS i F i 1 t_t 2< J ) ) >Peak _va J ue_f 2  THEN 
PeaF  _va 1  ue_f  2* ABS <  F 1 1 t_f  2( J  > ) 

END  IF 
J  *  J  *  ’ 

NE  X T  - 
NE  >  '  > 


.  i  Env 


_f :  f  E  nve  ioc  detection  over  fract- 
1  i  c  r‘  a  I  interval  of  F"  1  for  f  2 
1 usinc  1/2  cycle  difference 

D i f  f  _p 1 us_ f  2*0 

D  ,  *  *  IK;  nuc  f2-0 

F  OR  E  *  1  f  c  ■*  I  TO  5*  I  f  2* 1 

B  i  ♦ f  f  2-F i 1 t_f 2  <  E  » - F ilt_f2(E*3«If2» 
If  D i <  f  _  f  2 >D i f  *  _p 1 us_ f  2  THEN 

•  ;  .  ij«  _  t  i  *  v  1 1  f  __  f  2 

ENL  I f 

I  f  L  .  f  f  ♦2/[)jff_**i  nu?  f  2  THEN 
L  ,  •  *  r  ,  n  o  -  f  c  *  1 1  •  ♦  f  i 


•value  of  Envelop  at  P  interval  for  f2 
If  I  .  *  f  _jj  i  u%  _  f  2  >ABS  <  D  i  f  f  _  i  nus_  f  2  >  THEN 

t  ’  .  ■  :  •  .  •  D  .  *  •  _p  1  ui  4  < 

r  r 

f've. oi  *2*Diff  t*inue_f2 

in:  :• 

- 1  r.^e  .  op  4  2  ” ;  E  nve  1  op_  f  2 

,f.N-  •  ; ' . r.  1  04  Tetf  DOP^t F R  F f F E C T ' F OR  REP 

AND  fe.  1  if 

]  ■  HR  1  '•  -  e  ;  oi  4  ’  !  'ftb'  1  f  nveloi  f*1  THEN 


FOR  E*lf  1  +  1  TO  S* I f  1  -*■  1 

D2*Filt_f1<E*2)-Filt  f  1  (E  +  1 ) 

Dt  *F 1 1 1 _ f 1 (E*1  )-F  1 1 t_f 1(E) 

! PRINT  '*D2-~:D2.“DT-“:DT  ."(D2-D1  >*";D2-D1  ,"E-“;E 
IF  D2*0  OR  SGN(D1 )*SGN(D2>-0  THEN 
GOTO  S71 
END  IF 
NEXT  E 
Dsq*D2-D ! 

Ap*F i 1 t_f 1(E) 

Bp-D1 -Dsq/2 
Cp*Dsq/2 

Envelop_f 1 1»(Ap-(Bp‘2)/(4*Cp) ) 

•  PRINT  ‘Envelop_n 1-“;Envelop_f  11 

Ee*E-B p/(2*Cp) 

f PRINTER  IS  70S 

PRINT  "Ee_f1*";Ee 

PRINTER  IS  1 

Qc-0 

IF  Ee>4*I f 1 ♦ 1  THEN  Qc*-1 
IF  Ee<2*  I f 1 +1  THEN  Qc-  +  1 
ELSE 

FOR  E-IF2+1  TO  S*If2+1 

D2-F i 1 t_f2(E+2 ) -F i 1 t_f 2<E+1 ) 

DI'Filt  f2(E+1 )-Fi 1 t_f2(E ) 

(PRINT  *D2«“  :D2,  “D1  - “  :D1  ."(D2-D1  >*":D2-D1  ,  "E-";! 
IF  D2-0  OR  SGN(D1 )+SGN(D2)-0  THEN 
GOTO  626 
END  IF 
NEXT  E 
Dsq*D2-D1 
Ap*F i 1 t_f 2(E) 

Bp-D1 -Ds/2 
Cp-Dsq/2 

Envelop_f 22* ( Ap- ( Bp  *  2 ) / ( A*Cp ) ) 

(PRINT  "Enveiop_f22*";Enveiop_f22 
Ee-E -Bp/ ( 2*Cp ) 

(PRINTER  IS  70S 
PRINT  "Ee_f2-’;Ee 
PRINTER  IS  1 
Qc-0 

IF  Ee>4*If2+1  THEN  Qc*-1 
IF  Ee<2»If 2*1  THEN  Qc-*1 
END  IF 

IF  Xx-1  THEN 
(PRINTER  IS  70S 
Ror  *Ror  +  Qc 
PRINT  “Ror-":Ror, 

Rot r (P ) -Ror 
Qccr(P)*Qc 
PRINT  "Qc-";0c 
Qr-Qr+Qc+K  1 
PRINT  " Q r * " ; Q r 
PRINTER  IS  I 
ELSE 

Rob*Rob«-Qc 
PRINT  "Rob*" ;Rob 
Rotb(P)*Rob 
Occb ( P  > *Qc 
PRINT  "Qc*” ;Qc. 

Qb*Qb+Qe+K1 


*?> 


PRINT  "Ob* " : Ob 

END  IF 

E_F 1 *E nv* i op_F  1 
£_F2-4*Enve 1  op _F2 

PRINT  “Enve lop_f  2/Enve  J  op_F  1  * “ ; E nve lop_F  2/Envelop_F 1 

! . 

_t : 1 Compu t a t i on  oF  the  anole  From  the 
!  arc t an ( £ nve 1 op_f 2/E nve l op_F 1  J  ana 
Icounting  oF  the  whole  quadrants 

i . 

BEEP  SO 0  ,  .  2 .  . .  •  •  •  . . 

PRINTER  IS  1 
DEG 

IF  Xx*1  THEN  I  Arc  t  For  Red  laser 

IF  E_F 1-0  THEN  'Tate  care  oF  divide  by  tero 
IF  Qpsr-1  OR  Qpsr*3  THEN 
Tar- (Qsr )*90*90 

GOTO  716 

else 

Tar- (Qsr  > *90 
GOTO  716 
END  IF 
END  IF 

Ra*E_i 2/E_F 1 

IDetermine  Quadrant  &  1 o ta 1 _anale< P ) _ 

IF  ftESlRa ) >-60  THEN 
IF  SON i Ra ) • 1  THEN 
Ra-60 
ELSE 
Ra  -  -60 
END  IF 
END  IF 

IF  SGN(Ra>-1  THEN  IPlus  Ratio 
IF  SGN(E_F1)*1  THEN  •  Deno  plus 
Opt  r- 1 

else 

Opt  r-3 
END  IF 

Tar-(Qt (Qpsr.Qptr ) ♦Qsr )«90*ftTN(Ra ) -900 

ELSE 

IF  SGN(E_t 1 )* 1  THEN 
Qptr-4 
ELSE 
Qp  t  r  -  2 
END  IF 

Ta r  *  <  Ot  <  Qpsr , Qp t  r ) *Qs r ) *90*90*ATN( Ra ) -900 

END  IF 

Qsr*Ot(Qpsr,Qptr)+Qsr 
IF  P>1  THEN  Qpsr-Qptr 
PRINTER  IS  70S 

PRINT  P,PR0UND< Tar,-2) .Qpsr.Qptr.BT  <Q),B1  (P) 

PRINTER  IS  1 
T a_red ( P  > * Tar 
PRINTER  IS  70S 

PRINT  P,"Ia_red-";PR0UND(Ta_red(P).-2) ,’E_Fl-’;E_f1 .•E_E2-";E_F. 
PRINTER  IS  1 

ELSE!Arc_t  For  blue  laser 

IF  E_f1*G  THEN  ITare  care  oF  divide  by  tero 
IF  Opsb* 1  OR  0psb*3  THEN 
Tab- (Osb ) ♦90 ♦9u 
GOTO  763 


ELSE 

T  ab* ( Qsb ) *90 
GOTO  763 
END  IF 
END  IF 

Ra-E_f2/E_f1 

(Determine  Quadrant  A  To t a  1 _anq le ( P )_ 

IF  ABS l Ra ) > *60  THEN 
IF  SGN( Ra ) *  1  THEN 
Ra*6C 
ELSE 
Ra* -60 
END  IF 
END  IF 

IF  SGN(Ra)-1  THEN  'Pius  Ratio 
IF  SGN(E_f1)“1  THEN  !  Deno  plus 
Qptb* 1 
ELSE 
Qp  tb- 3 
END  IF 

Tab* l Qt  t  Qpsb , Qp  tb ) +056 )*90*ATN(Ra) -900 
ELSE 

IF  SGN ( £_f 1 ) -  1  THEN 
Qptb*4 
ELSE 

Qptb* 2 
END  IF 

Tab* ( Q t  <  Qpsb , Qp  t b ) +  Qsb )«90*90*ATN(Ra>-900 
END  IF 

Qsb'Qt ( Qpsb .Qptb ) +0sb 
IF  P>1  THEN  Qpsb*Qptb 
PRINTER  IS  70S 

PRINT  P,PR0UNDUab,-2>. Qpsb. Qptb. BUQl.BT  tP> 

Ta_b 1 <  P ) *  Tab 

PRINT  P . ■ Ta_b 1  -  * ;PR0UND< Ta_b 1 ( P ) , -2  > , "E_F 1  *  * ;E_f 1 . *E_F2- ' ;E_f 2 
END  IF (For  red  or  blue 
PRINTER  IS  1 
NEXT  X* 

RAD 

PRINT 

(This  completes  the  angle  computation 
(For  the  red  and  blue  for  this  K1  interval  P 
P-P+  1 

(Test  to  increment  Qo  or  not 
If  P-G  THEN 

Qq*0 

else 

Rmdr*P  MOD  199  !P  starts  at  U 
IF  Rmdr*0  THEN  Qq-Qq+1 
END  IF 
(NEXT  Q 

IF  Q > • Nnn- N - 2S - S t a r t  THEN  RETURN 
GOTO  363 

Disp_out : (Graph  data 

Dpcr*.064  (Degrees/count  gimbal  red 

Dpcb*.064  (Degrees/count  oimbal  blue 

DUMP  DEVICE  IS  70S 
PRINTER  IS  70S 
PEN  2 

PRINT  CHRS(12';  (Clear  alpha  A  form  teed 
DISP 


VjV.V 


**7** 


■'% 

797 

7Sb 

79s 

80U 

801 

802 

802 

804 

809 
80b 

80  7 
808 
80  4 

810 

81  1 
8 1 : 

f  1 

t-H 

819 

8 1  b 
817 
918 
■■  i  9 

82  0 
821 
82.' 
824 
H  jr 

82b 

829 

82i 

831 
833 
8  34 

839 
83b 
837 
83* 
8  33 

840 
34  1 

9. 

9 

847. 
<2 
^  u> 
u  ' 
84- 

94'' 


b  9 1 
89  - 

'-.U 
c.  r 

89b 

U 

-'l  E 


GIN  IT  1  I n 1 1 1  a J 1 / e 

GRAPHICS  ON  IRaster  or 

FOR  Oqq‘0  T 0  Qq  'Graph  j'de* 

FOR  R*  1  TO  2  !red*1.  blue*2 
PEN  2 

WINDOW  -2S , 200 . -S00 ,SOO ! S«t  Hindou 
AXES  20.100 
LORG  6 

FOR  1-0  TO  180  STEP  40 
HOVE  1.0 
LABEL  0qq*200 4  1 
NEXT  1 

I 

LORG  S 

FOR  I* -400  TO  400  STEP  100 
HOVE  U.I 
LABEt  I 

NEXT  J 

LORG  9 

MOVE  100.479 
PEN  b 

LABEt  "ANGLE  AS  A  FUNCTION  Of  ENVELOP  SAMPLES" 
HOVE  180.429 

ir  r *  1  then 

PEN  £ 

LABti  "RED' 

HOVE  180.  400 
LAEiEl  NdfieS 
PEN  b 
END  If 

if  r*2  then  label  "Blue  * 

HOVc  180.  40 u 
LABEL  NameS 
PEN  2 

HOVE  -18.0 
lQRO  b 
LABE  L  "Samp' 

LORG  9 

HOVE  IDT1.  -  4  7l 

lALEl  T it  ieS; “ENVELOP  SAMPLE  RATE.  3KH/* 

PEN  0 
HOVE  0 . 0 
EE  N  2 

f  CjR  X*i  1 0  199  'urapo  cata 
IE  R  *  1  THEN 

T!  "  1  red  uq  ,  4  20N  4  X  4  4  1  T  j  ^ 0  <  4  1 
L 1  3E 

Tb-Ta  M  <Qqa«200*X«4  1  Ta  bi'4> 

ENG  I* 
fiRAW  X  .  Tt 
NE  >  > 

Pt  N  0 
MOVE  0.0 
P{  I.  1 

E OR  X-0  T [ 1  199  'Graph  a jmha I 
I*  R  •  THE  t. 

Bbb- B 1  <  Qqq»2004  X  4 4  )  »Dpr  r  B 1  '  4  1  «Dpc  r 
Bt’b'  •  E  2  '  (jqq*?0li  ♦  X  4  4  t  *Uprb  82 ' 4  1  *Dpcb 
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9b  b 

DRAW 

X .  BbU 

86  i 

NEXT  X 

86. 

PEN  0 

dto 

PAUSE 

864 

! DUMP  GRAPHICS 

86b 

! PRINT 

CHRS (12) 

866 

GCLEAR 

86. 7 

NEXT  R 

86o 

NEXT  Qqq 

869 

ALPHA  ON 

870 

PRINTER  IS 

1 

871 

DISP  "STOP- 

6  7. 

PAUSE 

673 

RE  TURN 

874  Checkstat:  'Status  check  on  Buffer. 

8>S  PRINT  -  REGISTER  STATUS  #Buffer" 

876  FOR  X*i  TO  6 

67/  STATUS  eBuf  f er , X ; A 

6’b  PRINT  ‘REGISTER  ‘;X.A 

878  NEXT  X 

y*v  'DISP  ’PAUSE  PRESS  CONTINUE  TO  CONTINUE" 

o n  Bt  t y 

66.  'PAUSE 

6b  HP  :  ’  1 

66a  PRINT  LHR$ ( 1 > 

86'  RE  TURN 

6fc’t  t  r  r_  r  out  :  'Routine  tor  error  recove  r  y 

66  7  BEEP 

860  BEE1 
BEE' 

881  DISP  ERRMS 

861  PAUSE 

68.  GOoUE  Disi-'_out 
68  PAUSE 

894  ! 

8%  Pile  Oats!  'Transfer  of  data  from  Buffer 
886'  'to  dist 

88’  ! 

898  BEEP 

89A  DISP  ‘INSERT  a  FORMATTED  DATA  DISK! !'•'!'  -  PRESS  CONT  TO  CONT‘ 

9'  ■  Pa j _,t 

901  ! C reate  a  BDAT  file  for  the  data. 

SOU  INPUT  ‘NAME  Of  DATA  ERE  TO  BE  CREATED7"  .Name* 

8i"  IE  Na*el-“  THEN  9d. 

-<  0  Bt  if 

90S  DISP  ‘CREATING  A  DATA  FILE" 

Si.’ •  !  Record  lenatf  i<  1 

907  CREATE  BOAT  NaeeSA  * : INTERNAL ",8«Nnn*4.2 

908  lAsstgn  an  I/O  path 

909  ASSIGN  #f i le  TO  NaeeSA" : INTERNAL " 

RIO  GOSUB  Check_stat_r 

91  1  PRINT  CARS'  1? i 

-•1.  DISP  "  TRANSf  ERRING  DATA  FROM  BUFFER  10  F  J L  E  " 

9 1  -t  OUTPUT  *F  I  ie'.Br  I  «  )  ,B1 »  •  >  .B1<  •  )  .B2<  •  t 

814  ASSIGN  %F i le  TO  • 

9 1 S  HA  I T  S 

-IE  DISF  "TRANSFER  COMPLETE" 

17  HA  I T  S 

-’6  DISP 

'  ■-  PRINT  CHRU  12  ' 

<.(  RETURN 


69 


R 


921 

922 

923 

924 

925 

926 

929 

930 

931 

933 

934 

938 

939 

940 

941 

942 

943 

944 
949 
946 
94  7 

948 

949 

990 

991 

992 

993 

994 
999 
996 


File_array1  ITransfer  fro*  file  to  array 

; 

INPUT  “NAME  OF  DATA  FILE7" .NameS 
ASSIGN  §F i ie  TO  Na*e$4“ : INTERNAL “ 

BEEP 

DISP  *  transferring  from  data  file 

ENTER  #File;Br<»),Bl<*),B1iO,B2 IM¬ 
PRINT  CHRSt 12) 

DISP  ‘TRANSFER  COMPLETE*' 

WAIT  3 

Check _s tat 
Check_stat_r 
Check  stat 


TO  AN  ARRAY* 


IGOSUB 
IGOSUB 
IGOSUB 
DISP 
RETURN 

Check_stat_r : 

BEEP 
PRINT  “ 

FOR  R-1  TO  9 

STATUS  «<File.R:A 
PRINT  ‘REGISTER 


(Status  check  on  File. 

REGISTER  STATUS  #File“ 


NEXT  R 

BEEP 

WAIT  9 
PRINT 
RETURN 
END 


‘  ;R  ,A 

!  DISP  “PAUSE  -  PRESS  CONTINUE  TO  CONTINUE* 


1  PAUSE 
CHRS< 12) 


Cross  Reference 


<  e  e  \ 


»  Numeric  Variables 


H 

877 

878 

34/ 

948 

8_to_a_rese  t 

Ha 

306 

307 

312 

320 

321 

326 

Hr. 

278 

281 

287 

He 

28 1 

28  3 

Hr 

672 

6/7 

627 

630 

B  ' 

2 1  0 

266 

381 

412 

719 

766 

866 

913 

9  30 

t i _peat 

"  ■  K  ' 

363 

381 

412 

Bw 

21  1 

266 

382 

41  3 

868 

913 

9  30 

62  peak 

6  'k  1 

370 

382 

413 

Bo 

307 

3 1 2 

321 

326 

Bt.2 

Bbb 

866 

368 

860 

B  1 

208 

268 

386 

498 

913 

930 

Bt 

6  78 

6  7/ 

680 

628 

630 

63. 

f  r 

208 

264 

387 

490 

913 

930 

bu*ter  array 

Pu  1  dDe  i 

383 

38  7 

396 

410 

486 

490 

498 

613 

'  neck  si,)' 
Pec*  vlat  r 
l  omp  a  no  .  ►- 

or  t  . 


1,  ?[ 

C  7  7 

680 

629 

6  30 

832 

-  r  •. 

r  /  * 

bbl 

Sb  ' 
66  3 

671 

671 

6/4 

6  1  b 

6  ’  9 
621 

62' 
62  b 

•  •  ♦ 

♦  ' 

426 

427 

428 

430 

431 

i.7  i  ♦  ♦ 

♦  , 

629 

630 

631 

633 

634 

l)|M 

m i nus  ♦ 1 

424 

430 

431 

438 

441 

r.  .  , 

r  i  nu  ♦  . 

62  7 

6  i  i 

694 

641 

644 

.  i  *  ♦ 

f  l  us  t  ■ 
p  l  U  V  ♦  . 

42  < 

42  / 

428 

438 

439 

0  i  ♦  ♦ 

628 

630 

631 

641 

642 

I  P  /  on 


.  ■  t  Oi1  i 
I  '  I  U<  t  a 
iSP  <1a  1  a  . 

or.  i  t  ran  ■•*  ► 

.  ,  .  ’  a i i  HSH 

l,‘p  /8-  8S». 

,  K’H 


:  *  • 

'  7  4 

7  1 

t.2k 

P  "« 

o.  i 

<4 

C 

%  '  *9 

l,  IS 

668 

6b  it 

1  . 

ct  /  • 

6HI 

b  ;  / 

bib 

b  1  lJ 

b.H 

62/ 

632 

bKP 

►.m 

/II  1 

/or 

7  >( 

729 

7  3b 

SHI- 

k  .  1 

l.Hl 

( 

6  HO 

’  K  ■ 

t.  <. 

P  P 

»  '  - 

b  4 1 

*  .  ■  'i  *  '  4  "4 

r 

4  4  • 

i  : 

r  M 

P  'I 

»  .*1  * 

-.Ml-  ♦  .  l«'». 

1 . 40 

i  t 

U.  < 

►  /(I 

71 


273  281 


Epa 

Err.rout 
F  1 
F2 
Fb 

F i le.ar  ray 
F i le.aata 
F i 1 t.bw.i nt 
Frit.fl 
F  j lt.f 1  lb 
r 1 1 t.f 1  1  r 
F  i  1 1  _  ♦  1  b 
F i i t.f 1 r 
F i 1 t.f  2 
Fi i t.f 22b 
r 1 1 t.f  22 r 
F i 1 t.f  2b 
F  j  ]  t  _  f  2  r 
Gal  cyc.frg 
Of 

Op  i  o 
Hp  IP 

I 

T  .  ' 


i  ♦  / 


I 

•  I 

y 

K a > _ >  axe: 

Inn 

hr 

Nul t.f un 
N 

Nnn 

Us 

9un  sarip 
9um_samp let 

Ot 


F«a»  _va 1 ue. da  t  a 
~>ak  _va  1  ue.  ft 
: ea^  value  (  '  1 1 
•  eax  _  vd i ue  *  '  1  r 
F*ax  _va 1 ue_f 1 1 
p«ak  _va 1 u*_f  1  r 
c>a6  _va  1  ue_  *  2 
F’eak  _ wa  1  ue_f  22t 
F'eatr  i ue_<  22  r 

F*ak  _va 1 ue  f  2t 
F#alr  _va  1  ue_  *  2  r 

:  r oduc  t  f  1 


297 

<  -DEF 

303 

312 

387 

395 

297 

<  -DEF 

317 

32b 

490 

498 

250 

251 

307 

321 

266 

402 

406 

407 

426 

560 

561 

572 

266 

506 

509 

510 

529 

618 

619 

627 

277 

281 

87 

281 

307 

321 

83 

196 

197 

198 

201 

202 

204 

220 

221 

84 

804 

805 

806 

807 

810 

81  1 

812 

813 

262 

266 

305 

306 

38  3 

367 

395 

425 

426 

c  c 

594 

595 

264 

26b 

319 

320 

48b 

490 

498 

528 

529 

i  1  7 

63Q 

640 

378 

402 

406 

407 

409 

483 

505 

509 

510 

51. 

88 

186 

266 

346 

380 

426 

485 

651 

660 

89 

9f 

19  '-:  • 

■  -DEF 

1  94 

206 

204 

210 

2M 

386 

387 

392 

394 

395 

400 

489 

490 

495 

49' 

«96 

50  3 

279 

283 

285 

26' 

305 

306 

319 

320 

348 

386 

394 

489 

44' 

7*5 

215 

253 

254 

255 

256 

785 

907 

1bc. 

167 

1  bG 

25  j 

279 

280 

i  [ 

~>0( 

312 

3 i  4 

3  1  9 

320 

32b 

326 

'<  a 

6  ' 

4  1  3 

64/ 

64b 

65  7 

66b 

71b 

719 

72. 

7.6 

766 

767 

766 

776 

778 

781 

2  76 

379 

406 

407 

AM 

509 

510 

:6  ' 

986 

94b 

72 


Product  * 2 

490 

491 

498 

499 

Q 

365 

367 

380 

48S 

719 

766 

78S 

Gt 

363 

367 

660 

661 

k 

c,  93 

S94 

59S 

638 

639 

640 

64S 

648 

649 

6S0 

BBS 

658 

6S9 

660 

Qccb 

2S6 

6S8 

Gcc  r 

2S6 

648 

Ci  j 

Clpsb 

730 

7S3 

760 

762 

763 

766 

Upst 

683 

706 

713 

7 1 S 

716 

719 

Gptb 

749 

75 1 

753 

7S6 

7S8 

760 

762 

763 

766 

Gp  t  r 

702 

704 

706 

709 

711 

713 

7 1 S 

716 

719 

k 

3S6 

779 

782 

798 

Uqq 

798 

806 

84S 

847 

8S6 

8S8 

868 

Jr 

361 

362 

36  5 

6S0 

6S1 

Gsb 

144 

731 

734 

7S3 

760 

762 

Gsr 

143 

684 

687 

706 

713 

7 1 S 

Gt 

1  34 

<  -DE6 

1  35 

136 

137 

138 

139 

140 

141 

u; 

706 

713 

7  IS 

7S3 

760 

762 

Uuo  t ier  t 

34c. 

346 

\ 

7qc, 

820 

829 

844 

8S5 

867 

946 

947 

948 

946 

Pd 

691 

693 

694 

69S 

697 

700 

706 

713 

73b 

741 

741 

742 

744 

747 

7S3  _ 

760 

Ranoestored 

Tia  r 

f  t 

r 

78'.: 

Pop 

6SC 

6S6 

6S  i7 

lor 

645 

646 

647 

b  » 

186 

Pott 

266 

667 

P  o  t  r 

2  St 

647 

•9, 

208 

209 

210 

21  1 

213 

21S 

• 

26  ' 

28* 

--1  S  P 

tart 

14  M 

355 

361 

786 

Suri_i  1 

384 

388 

396 

402 

U'r-  ' 

qq  i 

499 

SOS 

>y  df-  i  _t  . 

nys_.de i  _  red 

Ta_b  1 

188 

76  7 

768 

847 

Ia_  rec 

1  RP 

722 

7  2C 

84‘ 

"  at 

:  1 

7  14 

753 

760 

766 

767 

■  a  r 

6P4 

687 

706 

713 

719 

722 

't 

p4< 

847 

849 

e  -  t  a  a  t  i 

*  c  Pt  X 

’ot. anaie.pea* 

>4  1 

_  ii 

208 

2U9 

210 

211 

214 

280 

28. 

28b 

3b 

391 

382 

386 

38  7 

394 

395 

41 1 

485 

489 

4hi, 

497 

498 

6U 

843 

846 

847 

849 

8S0 

8S4 

85b 

86c 

86  r 

861 

876 

877 

87p 

879 

* 

IV  II 

• 

iv\ 

7PC 

4 ' 2 

413 

488 

64. 

t  ' 

* 

28. 

2t_- 

286 

28t 

*  v  h  r  i  <9  C  •  .  f ' 

-  1  : 

'  t. 

1  V 

•  <*  "i*1  i 

-  • 

% 

M  i 

Ml  ' 

M(!>h 

M21 

r  ,  f  : 
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-•  I/O  Path  Names 

*BuF 

«  B  u  f  f  e  r 

s»F  1  le 

SGpio 

•  Line  Labels 
A_to_d_rese  t 
8rc_t 

BuF  f er_ar  ray 
Bw_labe 1 
Check_stat 
Check  _sta  t_r 
Comments 
Comp_ang le 
Con_t 1 
Co r._f  ' 

Constants 
Disp_con 
D  i  sr_data 
D  isp_da  t  a_1 
Disp_data_2 
Disp_out 
Or,a_t  ranst  e  r 
Dopp ler_cor 
E  nv_t  1 
E n\ _f 2 
i r  r_rou  t 
f i ie_ar  ray 
7 i le_data 

r  irt 

1  nde  « 

Initialize 
ha  i r 

1u  1  t_ ♦ un 
hur  samt 
Hanae_s  to  red 
Test 

Te^  t  da  t  a 

va  r_  i  1 1- 1 


194 

200 

285 

877 

909 

913 

145 

206 

672 

<  -DEF 

290 

874 

910 

943 

19 

<  -DEF 

176 

336 

374 

<-DEF 

479 

<  -DEF 

81 

179 

159 

190 

420 

<  -  DEF 

52 1 

<  -DEr 

11 

886 

167 

9  2  c 

160 

895 

1 55 

246 

74 

'  -  DEF 

148 

<-  -DEr 

'  74 

29.' 

150 

182 

268 

<  -  DEF 

206 

914  925 


<-DEF 

<-DEF 

<-DEF 


788  <-DEF 
<  -DEF 


<  -  DEr 

<  -DEr 

<  -  DEF 

<  -  DEF 


<  -DEF 
^  -DEr 


L i ne  Numbe  r 


•PROGRAM  FOR  DA T A  ACQUISITION  AND 
1  DA  r  A  PROCESSING  or  RANGE  DATA. 


1  ^  i  c  prooran  uses  the  difference 
•method  tor  detemninQ  the  amplitude 
'•or  ♦ '  A  * 2 . 

1  0  r  apr.  i  no  is  included  for  200  enve- 
1  1  opes / g r apn .  *  ot  graphs  depends  on 
•the  number  os  Num  samples. 

•arameteri  o*  consideration:  Ratios 
•c*  fD'*'.  degrees/count  red(Dpcr), 
iThe  phase  correction,  doppler 
1 e  *  *  ec  t  is  inc  1 udec . 

'The  parabolic  fit  is  used  to  deter- 
•mine  the  location  of  the  peat  using 
1  *  *  *  area^er  o’  the  amplitudes  ot  f’ 

'6  ♦.  .  Tue  phase  correction  is 
.noepenoe’-*  »or  the  red  A  blue.  — 

'AD  Interface  Switch  position  NORMAi 

r  r  I N  *  t  R  I '  * 

PR  IN'  CAPS 

JR  IN'  •'AS*'  '  S ,  S  .  ‘RANGE  PROGRAM' 


IIS  GRAPHICS  OFF 
t16  RAD 

117  ON  KEY  4  LABEl  ‘GRAPh". GGSub  Diw  o-.' 

118  ENABLE 

119  Gp 10* 1 2 

120  Hpib*704 

121 

123  Constant s  :  ' 

124  ! _ _ 

1 26  SF  *36000  ! Sample  Frequency 

127  GF  • 3000  !ga  1 vo  F requency  set  by  clock 

126  K1*12  (samples  at  <  1 

129  K2*6  (samples  at  F? 

130  K3*4  (samples  at  F3 

1  3 1  Sys_del_red* -0 ( -3  For  data  acquisition 

132  ‘*0  For  Test_data 

134  Sys_de )_b 1  *  -  0  !-3  For  data  acquisitior 

1 35  1 ♦ 0  For  T*s  t  _da  t  a 

137  (Sianal  aenerator  t  repuenc,  S  76 . 000*  i 

136  T I::::::::::::::::::::::::::::::: : 

139  Va r_ 1 i s t : ! Var lab le  list 

140  ! _ 

141  !Num_samp  :  number  oF  samples  to  ta*e  at  36*;  rata 

142  !Fb  1  Filter  bandwidtt 

143  ! N  :  INT  <  SF /Fb 1 

144  !Ns  ; INT  (Num_samp les/4 > 

145  !Br(  1  :  array,  saapled  data,  red 

146  ! B i (  l  :  array,  saapled  data,  blue 

147  !IF1  :  interpolation  Factor  For  FI 

148  *  I F 2  :  interpolat ion  Factor  For  F2 

149  ! F 1 (  )  :  convolution  muitipling  Function  For  FI 

150  !F2(  1  :  convolution  aulttpiing  Function  Fcr  <. 

151  'Start  1  starting  point,  deterained  by  N 

152  !0r  :  starting  saapled  data  point  For  next  K  1 

1  S3  (interval,  including  doppler  correction,  red. 

154  !Qb  :  same  as  Or  except,  blue 

155  !Xx  :  Xx-1,  loop  For  red;Xx*2  loop  For  blue 

156  ! B 1 t  )  1  spare  input 

157  !B2<  )  .*  spare  input 

158  !Filt_F1(J)  :  convolution  eleaents  in  K1  interval.  FI 

159  !Filt_F2(J)  1  convolution  eleaents  in  K1  interval.  F2 

160  ! B 1 ( P )  :  average  oF  Bit  >  in  K!  interval 

161  !B2(P)  :  average  oF  B2(  )  in  k2  interval 

162  !P  :  K1  interval  oF  coaputation 

163  !Ror  :  running  sum  oF  doppler  correction,  red 

164  !Rob  :  running  sum  oF  doppier  correction,  blue 

165  !Rotr(P)  :  running  sun  stored  as  a  Function  oF  P 

166  (Rotb(P)  :  running  sun  stored  as  a  Function  oF  P 

167  !0ccr(P)  :  doppler  correction  For  respective  K1  interval, 

168  !0ccb(P)  :  doppler  correction  For  respective  K1  interval, 

169  !0psr  J  quadrant  past,  red 

170  (Qpsb  :  quadrant  past,  blue 

171  IGptr  :  quadrant  present,  red 

172  IQptb  :  quadrant  present,  blue 

173  !0sr  :  running  sum  oF  quadrants,  red 

174  !0sb  :  running  sum  oF  quadrants,  blue 

175  !Tar  :  total  angle  red,  K1  interval 

176  !Tab  :  total  angle  blue,  K1  interval 

177  !Ta_red(P)  :  total  angle  red.  P  K1  interval 

178  !Ta_bl(P)  :  total  angle  blue.  P  K1  interval 

179  (Whole  quadrant  update  table 


red 

blue 


•l 

i 

A 


3 


■J 


D  I  H  U  t  '  l  ,  4 

(JtM.il-' 

(Jt  It  .  M  ,  •  ’ 

0 1  <  j  .  1  •  •  ' 

G  t  (  2 .  3  >  •  1 
Qt  <  3.2  - 1 

Gt i 3 . 4  i  *  1 

(jt  i  4  .  '  i  * ' 

Gt 1 4  .  j  I  *  1 

Gsr-10  'Quadrant  sun  red  starts  at  10 
Gsb*10  (Quadrant  sue  blue  starts  at  10 

ASSIGN  frGp i o  T 0  12 

!PRINT  “  REGISTER  STATUS  BEFORE  DMA' 
1 GOSUb  Chec*i_stat 

I  ; 

Main:!  Main  program 


PR  It/  CHRS  M  2"  i 
GOSUb  A_tc_d_ reset 
GOSUb  Range_stored 
6EEP 

GOSUB  Mu  It  fun 


(Generate  Multiplier 
(functions. 


4 


GOSUb  Co«p_ana]e  (Convolution,  envelop  detection,  doppler  correction, 

(A  Arc  tar  tor  red  and  blue 

DISP 

GOSUb  Disp_out 
GOTO  1 


Range  stored: ! 


s 

% 

I 

ft 


$ 

•2 

s 


INPUT  ‘RANGE  DATA,  Type  R/r  :  STORED  DATA'/,  Type  S/s",AaS 
IP  AaS-‘R“  OR  Aa$-‘r*  THEN 
GOSUB  Num_samp 
(GOSUB  Dma_transfer 
GOSUB  Test_data 
GOSUB  File_data 
GOSUB  Buffer_array 
GOSUB  Filt_bu_int 
RETURN 
END  IF 

IF  Aa$-‘S“  OR  Aa$-*s“  THEN 
GOSUB  F i le_ar ray 
GOSUB  Buffer_array 
GOSUB  Filt_bu_int 
RETURN 
END  IF 
BEEP 
BEEP 
BEEP 

GOTO  213 

Num_samp:  ! Input  the  total  number  of  samples 

! _ 

!  to  be  taken.  Maximum  of  32K 


INPUT  "TOTAL  NUMBER  OF  SAMPLES  TO  BE  TAKEN?" ,Num_samples 


243  ALLOCATE  INTEGER  Rotr(0:Ns> ,Rotb(0:Ns) ,Qccr(0.*Ns) ,Gccb(0:Ns) 

245  ALLOCATE  REAL  Ta_red( INT(Num_sanples/K1 ) ) , Ta_bl ( INT(Nun_sanples/K1 ) ) ,Rot( I 
NT (Num_samples/K1 ) ) 

246  {Create  an  un-named  buffer  so  greater  than 

247  (32k  can  be  collected. 

248  ASSIGN  ©Buffer  TO  BUFFER  (Nun  samples] 

249  RETURN 

250  Dma_transfer : ! 

251  CONTROL  ©Buffer. 4;0 

252  CONTROL  Gpio,0;2 

253  CONTROL  Gpio,3;0 

254  (SET  TRANSFER  LOW 

255  CONTROL  Gpio,2;3  (SET  CLTO  TO  0,  XFER  0 

256  BEEP 

257  DISP  *  WAITING  FOR  RANGE  TRIGGER" 

258  TRANSFER  ©Gpio  TO  ©Buf fer jCOUNT  Num_samples 

259  CONTROL  Gpio ,2;2  (SET  CLTO  A  CLT1  TO  1 

260  CONTROL  Gpio,1;3 

261  (PRINT  "  REGISTER  STATUS  AFTER  DMA' 

262  (GOSUB  Check_stat 

263  BEEP 

264  RETURN 

265  Fi 1 tbw_i nt :! Input  filter  banduidth  in  hertz 

266  BEtP 

267  INPUT  "WHAT  IS  THE  DESIRED  FILTER  BANDWIDTH  IN  HERTZ?', Fb 

268  N-INT(Sf/Fb) 

269  PRINT  CHRS< 1 2 ) 

270  BEEP 

271  INPUT  "Real  samples  f1“12,  interpolation  factor?", Ifl 

272  BEEP 

273  INPUT  "Real  samples  f2-6.  interpolation  factor?", If2 

274  BEEP 

275  ALLOCATE  REAL  Filt  fid :K1*If 1 )  ,Filt  f2( 1 :K1*If 2> 

276  RETURN 

277  Test:!  Generate  data  -  needs  to  be  modified  for  this 

278  •  (program 

279  BEEP 

280  INPUT  "What  is  the  Envelop  Peak  Ampl i tude?* ,Epa 

281  BEEP 

282  DISP  "  GENERATING  DATA" 

283  Peak_value_data'0 

284  Gal_cyc_f rg-200 

285  Ab-0 

286  ALLOCATE  INTEGER  Mmm< 0 :Num_sampies_ 1 ) 

287  FOR  X"0  TO  Num_samples- 1  STEP  4 

288  Ac-Epa*SIN( (PI* < +Ab) )/<6*Gal  cyc_f rg)  +  (2. 7)*SIN(2*PI*Ab*Gf /Sf ) ) 

289  FOR  Y-X  TO  X+3 

290  Mmml Y ) * INT ( Ac ) + 128 

291  (PRINT  Y,Mmm(Y> 

292  OUTPUT  ©Buffer  USING  B' ;Mmm( Y ) 

293  NEXT  Y 

294  Ab«Ab+ 1 

295  NEXT  X 

296  DISP 

298  WAIT  5 

299  RETURN 


300 

301  Mult_fun:  (Generation  of  the  multiplying  function. 

302  ! _ _ _ _ _ _ _ _ 

303  (Multiplyino  function  will be  l SIN  t  X ) / ( X )  )  *C0S l Y > 

304  (PRINTER  IS ' 705 


305  ! ALLOCATE  REAL  FI <0: (N*If 1 ) ) .REAL  F2(0: <N-If2) > 

306  DIM  FI ( 0:2000 ) ,F2( 0:2000 ) 

307  RAD 

308  BEEP 

309  IGenerate  1/2  multiplier  Function  for  FI 

310  BEEP 

311  DISP  "Building  multiplier  function  for  FI* 

312  Fl (0)-1 

313  ! PRINT  "FI (0)-":F1 (0) 

31 A  FOR  P-1  TO  N*If1 

315  Aa-(PI*P)/(N*If1) 

316  Bb-(Aa*2*Gf )/Fb 

317  t PRINT  "P«";P 

318  1  PRINT  "N-";N*If1 

319  ! PRINT  "(PI*P)/N-“:Aa 

320  (PRINT  "(Aa*2*GF)/FB-":Bb 

321  FI (P)-(SIN(Aa)/Aa)*COS(Bb) 

322  SPRINT  "FI  ( "  ;P; " ;FHP) 

323  NEXT  P 

324  DISP  “Building  multiplier  function  for  F2" 

325  IGenerate  1/2  of  multiplier  function  for  F2 

326  F2(0)»1 

327  ! PRINT  "F2f 0 ) - " :F2 ( 0 ) 

323  FOR  P-1  TO  N*If2 

329  Aa=(PI*P)/<N*If2> 

330  Bb= ( Aa*4*Gf ) /Fb 

331  IPRINT  "P-“ ;P 

332  IPRINT  "N-“:N 

333  IPRINT  "(PI*P)/N-“:Aa 

334  IPRINT  "(Aa*4*GF)/FB-":Bb 

335  F2(P)-(SIN(Aa)/Aa)*C0S(Bb) 

336  IPRINT  "F2( "  ;P;  *•)«•*  ;F2CP) 

337  NEXT  P 

338  RETURN 

339  I. 

340  I . 

341  IConvolution  and  envelop  detection 

342  !.at  K1  intervals  for  fl  and  f2. 

343  ! . 

344  I :::::::::::::::::::::::::::::::::::::::::::: : 

345  Comp_angle: IFi 1  ter ,  envelop  detection, 

348  Idoppler  correction  and 

349  larctan  for  red  6  blue 

350  I  _ 

352  Rad 

353  Tot_angle_peak-0 

354  P*0  !  Iru t lal  l  za ton  of  array  for  the  out- 

355  Iput  angle  for  red  and  blue  lasers. 

356  IDetermine  first  data  point  to  start  con- 

357  Ivolution  considering  the  value  of  N  and 

358  I  the  need  of  a  +  slope  for  fl  for  synchr- 

353  Ionization. 

360  Quotient- (N+ 1)  DIV  K1 

361  Start-(Quotient  +  1  )*12 

362  BEEP 

363  DISP  "  CRUNCHING  DATA" 

364  PRINT  "Start*" jStart 

365  Qq-0  (Initialize  graph  index.  Increment 
Gr-Star,^2iS$sia&l?_?89  computations 

371  Qb-Start*2+Sys_dei_bl 
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372  FOR  Xx-1  TO  2  ! 1-red, 2-blue 

374  IF  Xx-1  THEN 

375  Q-Qr 

376  ELSE 

377  Q-Qb 

378  END  IF 

379 

380  Con_fi:i 

381  {Convolution  over  XI  interval  for 

382  !  F 1 

383  ! _ 

384  J-1 

385  Peak_value_f 1 -0 

388  FOR  X-Q  TO  Q+K1-1 

392  FOR  C-0  TO  If  1 -1 

393  Sum_f1-0 

394  IF  Xx-1  THEN  ! Xx-1 ,Red:Xx-2,Blue 

395  FOR  M-X  TO  X+N 

396  Product_f 1 -Br  <M+ 1 )*F1 < (M-X)*If 1 -C+If 1 )+Br(2*X-M)*F1 ( <M-X)+If1+C) 

397  Sum_f 1 -Sum_f 1 +Product_f 1 

398  !  PRINT  Q; TAB< 1 0 ) :X ; TAB (20 ) ; C; TAB< 30 ) ;M; TAB(40 > ; J; TAB(45) ; ( (M-X)* 
Ifi+C) ; 

399  !  PRINT  TAB(50 ) ; ( <M-X)*If 1 -C+If  1 ) ; 

400  !  PRINT  TAB(60) ; (2+X-M-1 ) 

401  NEXT  M 

402  ELSE 

403  FOR  M-X  TO  X+N 

404  Product  f  1-BKM+1  )*F1  ( (h-X)«If  1-C+If  1  )  +B1  (2*X-M)*F1  ( <M-X)*If  1  +C> 

405  Sum_f 1 »5um_f 1 +Product_f 1 

406  !  PRINT  Q; TAB( 1 0 ) :X; TAB(20 ) ;C; TflB(30) ;M; TAB(40 ) ; J; TAB(45) ; ( (M-Xi 
+  1 1 1 +C ) ; 

407  !  PRINT  TAB(50) ; ( ( M-X )#If1 -C+If 1 ) J 

408  PRINT  TAB(60);<2*X-M-1> 

409  NEXT  M 

410  END  IF 

411  '  F i 1 t_f 1 ( J) -Sum_f 1 

413  IPRINTER  IS  705 

414  ! PRINT  "Filt_f1(";J;"-)‘;Filt_f1(J) 

417  IPRINTER  IS  1 

418  IF  ABS(Fi 1 t_f 1 ( J) ) >Peak_va 1 ue_f 1  THEN 

419  Peal _va 1 ue_f 1  - ABS ( F i 1 t_f 1 ( J ) ) 

420  END  IF 

421  J-J+1 

422  NEXT  C 

423  NEXT  X 

424  I:::::::::::::::::::::::::::::::::::::::::: : : : 

425  Env_f 1 : lEnvelop  detection  over  K1 

426  linterval  for  fl 

427  ! _ _ _ 

428  Dif f _plus_f 1 -0 

429  Dif f_minus_f 1 -0 

430  FOR  E- If  1  +  1  TO  5+If 1  +  1 

431  Dif f_f 1 -Fi 1 t_f 1 (E)-Filt_f1(E+INT(If 1+K1/2)) 

432  IF  Dif f_f 1  >Dif f_plus_f 1  THEN 

433  Diff  _plus_f 1 -Dif  f_f 1 

434  END  IF 

435  IF  Dif f_f 1 <Dif f_ninus  fl  THEN 

436  D l f f_mi nus_f 1 -D l f f_7l 

437  END  IF 

438  (PRINT  "e-":E,*Diff_f 1-":Diff_f 1 

439  SPRINT  "Dif f _pius_f 1 *" J D i f f_plus_f 1 
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440 

441 

442 

443 

444 

445 
448 

447 

448  ! : : 

449  Con 

450 

451  ! _ 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

♦If2+C) 

462 

463 

464 

465 

466 

467 

468 

469 

*If2+C) 

470 

471 

472  ... 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484  !  : : 

485  Env 

486 

487 

488  ! 

489 

490 

492 

493 

494 

495 

496 
«97 

498 

499 


! PRINT  *Dif 1  ■“ ;Dif f_einus_f  1 

NEXT  E 

(Value  of  Envelop  at  P  interval  for  fl 
I  Diff _plus_f 1 >ABS(Dif f_mnus_f 1 )  THEN 
i nve 1 op_f 1 * D i f  f _p 1 us_f 1 
ELSE 

Envelop_f  iDif  f_einus  fl 
END  IF 

f  2: ! . 

(Convolution  over  K1  interval  for  F2 

_ 

Peak  value  f2* 0 
FOR  K-Q  TO  Q+K2- 1 
FOR  C-0  TO  If 2- 1 
Sum_f 2*0 

IF  Xx*1  THEN  ! X x “ 1 ,Red:Xx*2,Blue 
FOR  M-X  TO  X+N 

Product  f2*Br(M+1 )*F2< (M-X)»If2-C:tIf2)*Brl2«X-M>#F2< (M-X)«If2*C> 
Sum_f 2*3um_f 2+Productjf  2 

!  PRINT  Q :  TAB ( 10)  ;X; TABOO )  ;C; TABOO)  ;M; TAB(40)  ;J; TAB(4S) ;  ( (M-X) 

!  PRINT  TABOO) ;  ( (M-X)*If2-C+If2) ; 

!  PRINT  TABOO)  ;<2*X-M-1 ) 

NEXT  M 
ELSE 

FOR  M-X  TO  X+N 

Product_f 2*B1 (M+1 )*F2( <M-X)*If2-C+If2)+Bl (2*X-M)*F2( (M-X )*If2*C ) 
Sum_f2-Sum_f2+Product_f2 

!  PRINT  Q;TAB(  10) ;x;TAB(20);C; TABOO) ;M;TAB(40);J:TAB(45);( (M-X) 

!  PRINT  TABOO)  ;UM-X)*If2-C+If2); 

!  PRINT  TABOO) ;  <2*X-M-1 ) 

NEXT  M 
END  IF 

F i 1 t_f 2 ( J ) -Sum_f 2 
(PRINTER  IS  705 

(PRINT  “F  i  1  t_f 2 C; J ;"*)’; F i  1  t_f 2 ( J ) 

(PRINTER  IS  1 

IF  ABS(F i 1 t_f 2 ( J) ) >Peak_value_f2  THEN 

Peak_value_f 2-ABS ( F i 1 1 _ f  2 ( J ) ) 

END  IF 
J-J+1 
NEXT  C 
NEXT  X 


_ f 2 : (Envelop  detection  over  fract¬ 
ional  interval  of  K1  for  f2 
(using  1/2  cycle  difference 

Dif  f_plus_f 2-0 

Dif  f_minus_f 2*0 

FOR  E*If2+1  TO  5*If2+1 

Dif f_f 2*F i 1 t_f2(E )-F l ] t_f2(EO*If  2) 
IF  Dif f_f2>Dif f_plus_f2  THEN 
Dif  f_plus_f2-Dif  f _f 2 
END  IF 

IF  Dif f_f2<Di f f_minus  (2  THEN 
Dif  f _mi nus_f 2*D  i  f  f_T2 
END  IF 
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SOD 

NEXT  E 

S0 1 

(Value  of  Envelop  at  P  interval  for  f2 

502 

IF  Ditf _p 1  us  f 2>ABS<  Di f f  «inus_f2)  THEN  i 

503 

Envelop  f 2- D i f f  plus  f2 

504 

ELSE 

505 

Envelop  f2*Diff  Minus  f2 

505 

END  IF 

507 

(PRINT  "Envelop  f2*:Envelop  f2 

510 

S11 

Doppler  cor : ! 

1  512 

i 

513 

i 

514 

!  CORRECTION  FOR  THE  DOPPLER  EFFECT  FOR  RED 

515 

!  AND  BLUE 

516 

1  . 

•  517 

IF  ABSlEnveiop  f 1 ) >ABS<Envelop  f 2 )  THEN 

518 

FOR  E * I f  1  ♦ 1  TO  5*If 1 ♦ 1 

518 

D2“F ilt  f 1 (E*2)-Filt  f 1 C E ♦ 1  ) 

520 

01  *F 1 1 1  f  1  (£♦!  )-Filt  f 1 (E ) 

521 

•PRINT  "D2-" ;D2. "DI :D1 , " (D2-D1 >•“ ;D2-D1 , “E - " ; E 

i  522 

IF  D2-0  OR  SGN(D1 )+SGN(D2>-0  THEN 

I  523 

GOTO  526 

!  524 

END  IF 

525 

NEXT  E 

526 

Dsq-D2-D1 

527 

Ap-Filt  f 1 ( E ) 

528 

Bp-DI -Dsq/2 

529 

Cp*Dsq/2 

!  530 

Envelop  f 1 1  * < Ap- ( Bp ‘2 ) /< 4*Cp) ) 

531 

•PRINT  "Envelop  f 1 1  - “ : Enve lop_f 1 1 

532 

Ee*E-Bp/(2*Cp) 

533 

(PRINTER  IS  705 

534 

! ! ! !  PRINT  "Ee  f1-";Ee 

»  535 

PRINTER  IS  1 

536 

Qc*0 

537 

IF  Ee>4*If 1 +1  THEN  0c--1 

538 

IF  Ee<2*If 1 -1  THEN  0c-+1 

539 

ELSE 

540 

FOR  E-IF2+1  TO  5*If2+1 

541 

D2-Filt  f 2( E+2  7  - F i 1 1  f2(E  +  1) 

1  542 

01-Filt  f2(E+1 )-Fi It  f 2(E ) 

543 

(PRINT  "D2-";D2,"D1-";D1 ,’(D2-D1 )-" ;D2-D1 . "E-B ;E 

544 

IF  D2-0  OR  SGNCD1 )+SGN(D2)*0  THEN 

545 

GOTO  548 

546 

END  IF 

547 

NEXT  E 

548 

Dsq'D2-D 1 

549 

Ap-Filt  f  2 ( E  > 

550 

Bp*D 1 -Ds/2 

551 

Cp»Dsq/2 

552 

Envelop  f22*(Ap-(Bp‘2)/(4*Cp) > 

553 

(PRINT  "Envelop  f 22-" ;Envelop_f22 

554 

Ee*E-Bp/(2*Cp> 

555 

(PRINTER  IS  705 

556 

!  ! ! !  PRINT  *Ee  f2-";Ee 

557 

PRINTER  IS  1 

558 

Qc*0 

559 

IF  Ee>4*If2+1  THEN  Qc--1 

560 

IF  Ee<2Mf2-1  THEN  Oc-+1 

561 

END  IF 

562 

IF  Xx  *  1  THEN 

i 
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563 

564 

565 
966 

567 

568 

569 

970 

971 

972 

973 
97^ 
979 

976 

977 

978 

979 

980 

981 

982 

983 

984 
989 

986 

987 

988 

989 

990 
591 

992 

993 

994 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

609 
606 
907 
608 
509 

610 
61  1 
612 

613 

614 

615 

616 

617 

618 

619 

620 
621 
622 
623 


Ror*Ror+Qc 

•Ml  PRINT  *Ror**  ;Ror , 
RotrlPl-Ror 
Gccr<P)-Qc 
IPRINTER  IS  705 
I !  PRINT  "Qc* " ;Qc 
PRINTER  IS  1 
Qr-Qr*Qc*K1 
I!  PRINT  "Of" ;  0  r 
ELSE 

Rob*Rob*Qc 
!  !  PRINT  •Rob-';Rob 
Rotb(P)“Rob 
Occb(P) *0c 
IPRINTER  IS  70S 
PRINT  *Qc-‘;0c 
PRINTER  IS  1 
Qb*Qb+Qc*K1 
PRINT  *Qb»* ;Qb 
PRINT 
END  IF 

E_f 1 -Envelop_f 1 
E_f2*Envelop_f2 

IPRINT  "Envelop_f 2/Envelop_f 1  *  ‘ 


•  I 

!  I 

!  !  ! 

MM 

ii  i  i 


!  !  ! 
i  1  1 

Arc 


; E  nve 1 op_f  2 /E  n ve 1 op_f 


_t : ! Computat ion  of  the  angle  from  the 
!arctan(Envelop_f2/Envelop_Fl )  and 
{counting  of  the  whole  quadrants 


BEEP  500,. 2 
PRINTER  IS  1 
DEG 

IF  Xx-1  THEN  I  Arc_t  for  Red  laser 

IF  E_f1«0  THEN  I  Take  care  of  divide  by  zero 
IF  Qpsr-1  OR  Qpsr-3  THEN 
Tar*  ( Qsr  M90+90 
GOTO  634 
ELSE 

Tar-(Qsr ) *90 
GOTO  634 
END  IF 
END  IF 
Ra-E_f 2/E_f 1 
IPRINTER  IS  705 
MM  PRINT  Ra 

PRINTER  IS  1 

IDetermine  Quadrant  4  Total_angle(P)_ 

IF  ABS(Ra) >"60  THEN 
IF  SGN(Ra)-1  THEN 
Ra-60 
ELSE 
Ra  * -60 
END  IF 
END  IF 

IF  SGN(Ra ) *  1  THEN  IPlus  Ratio 
IF  SGN(E_f 1 )“1  THEN  !  Deno  plus 
Qptr* 1 
ELSE 
Qptr* 3 
END  IF 

Tar* ( Qt ( Qpsr ,Qpt  r ) +  Gsr ) *90+ATN( Ra ) -900 
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1 

b34 

ELSl 

■ 

626 

If  SGN<E  f  1  )*  1  THEN 

H 

63p 

Qp  t  r  -4 

■ 

-..:• ; 

else 

R 

626 

Qptr*2 

w* 

629 

END  IF 

I* 

630 

Tar-lQt(Qpsr,Qptr)*Qsr)*90+90*ATN(Ra)-900 

631 

END  If 

K 

633 

Qsr-Qt (Qpsr ,Optr  )*Qsr 

633 

IF  P>1  THEN  Qpsr-Qptr 

■ 

634 

IRINTER  IS  706 

K 

636  !  !  ! 

!  PRINT  P ,PR0UND( Tar , -2 ) , Qpsr ,Qpt  r ,B1(Q),B1(P) 

K 

636 

PRINTER  IS  1 

K 

638 

PRINT 

K 

639 

Ta  red(P)-Tar 

K 

641  !  ! 

!  PRINT  P,  ”Ta  red-" ;PROUND( Ta  red(P) ,-2) , “E  F 1  - “ ; E  f 1 , "E_F2-" :E_f2 

■ 

642 

PRINTER  IS  1 

■ 

643 

ELSEIArc  t  for  blue  laser 

B 

644 

IF  E  fl-0  THEN  !  Tat  e  care  of  divide  by  lero 

c 

646 

IF  Qpsb-1  OR  Qpsb-3  THEN 

646 

T  ab- <Qsb >*90*90 

K 

64  7 

GOTO  679 

K 

646 

ELSE 

■ 

648 

T  ab- ( Qsb ) *90 

■ 

660 

GOTO  679 

R 

661 

END  IF 

662 

.  END  IF 

663 

Ra-E  f2/E  fl 

664 

IDetermine  Quadrant  A  Total  anale<P> 

666 

IF  A8S ( Ra ) >-60  THEN 

k 

666 

IF  SGN(Ra)-1  THEN 

I 

PC,  7 

Ra-60 

K 

668 

ELSE 

W) 

669 

Ra--60 

| 

660 

END  IF 

661 

END  IF 

1 

662 

IF  SGN ( Ra ) -  1  THEN  IPlus  Ratio 

663 

IF  SGN(E  f1)-1  THEN  !  Deno  plus 

664 

Qptb- 1 

666 

ELSE 

K 

666 

Qptb-3 

K 

667 

END  IF 

1 

668 

Tab- ( Qt ( Qpsb . Qptb)+Qsb)*90*ATN(Ra)-900 

r>  r. 

ELSE 

670 

IF  SGN(E  f  1)-1  THEN 

Bj 

671 

Qptb-4 

E 

672 

ELSE 

K* 

673 

Qptb-2 

IK 

674 

END  IF 

I 

676 

Tab- ( Qt ( Qpsb, Qptb ) +Qsb ) *90+90 +ATN( Ra ) -900 

676 

END  IF 

E 

m 

677 

Qsb-Qt ( Qpsb , Optb ) *Qsb 

678 

IF  P>1  THEN  Qpsb-Qptb 

R 

678 

! PRINTER  IS  706 

EE 

680 

PRINT 

Eh 

682  !  ! 

!  PRINT  P,PR0UND(Tab,-2),Qpsb.Qptb,BKQ),B1(P)  ' 

Ey 

683 

PRINTER  IS  1 

| 

686 

Ta  bl (P ) -Tab 

686  !  ! 

!  PRINT  P , " T a  bl * " ;PR0UND< Ta  bl(P),-2)."E  F 1  - “ ; E  f1,"E  F2-":E  f2 

E 

687 

END  IFIFor  red  or  blue 

K 

688 

PRINTER  IS  1 

f 

J 

1 

( 

_____ 

85 

| 

689  RAD 

690  NEXT  X* 

691  RAD 

692  PRINT 

693  (This  completes  the  angle  computation 

694  !for  the  red  and  blue  for  this  K1  interval  P 

695  P-P*1 

696  PRINT  P 

698  (Test  to  increment  Qq  or  not 

699  IF  P-0  THEN 

700  Qq- 0 

701  ELSE 

702  Rmdr-P  MOD  199  !P  starts  at  0 

703  IF  Rmdr-0  THEN  Qq-Qq+1 

704  END  IF 

705  !Q>Ns-48,  either  red/blue  will  exit  to  Main 

706  IF  Q>-Ns-48  THEN  RETURN 

707  GOTO  372 

708  Disp_out : (Graph  data 

709  DUMP  DEVICE  IS  70S 

710  PRINTER  IS  70S 

711  PEN  2 

712  ! (PRINT  CHRSC12);  (Clear  alpha  &  form  feed 

713  DISP 

714  GINIT  (Initialize 

715  GRAPHICS  ON  (Raster  on 

716  FOR  Qqq-0  TO  Qq  (Graph  index 

717  FOR  R-1  TO  2  !red-1,  blue-2 

7 1 8  PEN  2 

719  WINDOW  -7, SO, -500,600  (Set  uindou 

720  AXES  10,100 

721  LORG  6 

722  FOR  1-0  TO  50  STEP  10 

723  MOVE  1,0 

724..  LABEL  Qqq*200+I 

725  NEXT  I 

726  ! 

727  LORG  8 

728  FOR  I--400  TO  400  STEP  100 

729  MOVE  0,1 

730  LABEL  I 

731  NEXT  I 

732  ! 

733  LORG  S 

734  MOVE  25,475 

735  PEN  6 

736  LABEL  "ANGLE  AS  A  FUNCTION  OF  ENVELOP  SAMPLES" 

737  MOVE  18,425 

738  IF  R-1  THEN 

739  PEN  2 

740  LABEL  "RED" 

741  MOVE  18.-400 

742  LABEL  NameS 

743  PEN  6 

744  END  IF 

745  IF  R-2  THEN  LABEL  "BLUE" 

746  MOVE  18,-400 

747  LABEL  NameS 


86 


753  MOVE  25.-475 

754  LABEL  TitleS: “ENVELOP  SAMPLE  RATE.  3KHZ ' 

755  PEN  0 

756  MOVE  0 , 0 

757  PEN  2 

758  FOR  X*0  TO  50  Sgraph  data 

759  IF  R-1  THEN 

760  Tb*  T  a_redl Qqq*200*X+4 ) - Ta_red( 4 ) 

761  ELSE 

762  Tb*  T a_b 1 ( Qqq«200+X+4 ) - Ta_bl t  A ) 

763  END  IF 

7 6 A  DRAW  X ,  Tb 

765  NEXT  X 

766  PEN  0 

767  MOVE  0,0 

768  PEN  0 

769  PAUSE 

770  ! DUMP  GRAPHICS 

771  'PRINT  CHRS (12) 

772  GCLEAR 

773  NEXT  R 

77A  NEXT  Qqq 

775  ALPHA  ON 

776  PRINTER  IS  1 

777  DISP  "PAUSE" 

778  PAUSE 

779  RETURN 

780 

781  Err_rout:  'Routines  for  error  recovery 

782  ! _ 

783  BEEP 

78A  BEEP 

785  BEEP 

786  "  DISP  ERRMS 

787  PAUSE 

788  GOSUB  Disp_out 

789  PAUSE 

790 

791  Buf f er_array : ! Transf er  of  data  from  buffer  to  array. 

792  ! _ 

793  PRINT  CHR$(12) 

79A  DISP  "TRANSFERRING  DATA  FROM  BUFFER  TO  A  ARRAYS" 

795  CONTROL  ©Buffer.S;! 

796  FOR  X-1  TO  Ns- 1 

797  ENTER  ©Buffer  USING  "#.B":Zzz 

798  B r ( X - 1  )  -Z z  z- 1 28 

799  ENTER  ©Buffer  USING  "*.B";Zzz 

800  &1(X-1)-Zzz-128 

801  ENTER  ©Buffer  USING  *#,B’;Zzz 

902  B1 (X-1)-Zzz-128 

803  ENTER  ©Buffer  USING  "•,B";Zzz 

80A  B2(X-1 ) -Zz  z- 1 28 

805  SPRINT  X-1tBr(X-1),BKX-l),B1(X-U,B2(X-1) 

806  NEXT  X 

807  DISP  "TRANSFER  COMPLETE* 

808  WAIT  2 

809  PRINT  CHRS (12) 

810  RETURN 
911 


:$ 

ft 

■  a  I1 


,yj 


V* 


i'iV* 

»sv 

I 

V*V 


812  F 

813 
8M  ! 

815 

816 

817 

818 

819 

820 
821 
822 
823 
624 

825 

826 

827 

828 

829 

830 

831 

832 

833 

834  ! 

835  F 

836  ! 

837 

838 

839 

840 

841 

842 

843 

844 

845 

846  .  . 

847 

848 

849 

850 

851 

852 

853 

854 

855 

856 

857 

858 
959 
860 
861 
862  ! 

863  ft 

864  ! 

865 

866 

867 

868 
869 


ile_data:  (Transfer  of  data  fro*  Buffer 
!  to  disk 


DISF  ’INSERT  A  FORMATTED  DATA  DISK!  -  PRESS  CONT  TO  CONT' 
PAUSE 

INPUT  ’NAME  OF  DATA  FILE  TO  BE  CREATED?" , Na*eS 

IF  NaaeS* ’ *  THEN  818 

BEEP 

DISP  ’CREATING  A  DATA  FILE* 

(Record  length  is  1 

CREATE  BOAT  NaaeSA’ : INTERNAL ’ ,Nu»_saaplec ,  1 
(Assign  an  I/O  path 
ASSIGN  ©File  TO  NaaeSA" : INTERNAL " 

PRINT  CHR$< 12) 

DISP  ’TRANSFERRING  DATA  FROM  BUFFER  TO  FILE* 

TRANSFER  ©Buffer  TO  ©FilejCOUNT  Num  samples 
HAIT  5 

DISP  ’TRANSFER  COMPLETE’ 

HAIT  5 

DISP 

RETURN 


ile_array:  (Transfer  from  file  to  array 


BEEP 

PRINT  CHR$ (12) 

DISP  "INSERT  DATA  DISK  --  PRESS  CONT  TO  CONT* 


DISP  "INSERT  DATA  DISK  --  PRESS  CONT  TO  CONT" 

PAUSE 

CAT 

BEEP 

INPUT  ’WHAT  FILE?*. Name$ 

INPUT  "Rec/File  for  f i le?“ ,Num_samples 
GOSUB  240 

ASSIGN  ©File  TO  NaaeSA* : INTERNAL ’ 

(Reset  file  pointer  to  1 
CONTROL  ©F  i  1  e ,  5 ;  1 

(Reset  buffer  pointer  to  1 
CONTROL  ©Buffer,5;1 

DISP  "TRANSFERRING  FROM  DATA  FILE  TO  BUFFER  ARRAY* 
TRANSFER  ©File  TO  ©Buf f er ;C0UNT  Num_samples 
DISP  ’TRANSFER  COMPLETED" 

(Reset  buffer  pointer  to  1 
CONTROL  ©Buffer, 5; 1 
(Close  ©file 
ASSIGN  ©File  TO  * 

HAIT  2 

PRINT  CHR$( 1 2 ) 

DISP 

RETURN 


to_d  reset : 


(RESET  DATA  TRANSFER  ELECTRONICS 
CONTROL  Gpio,2;2  (SET  CTLO  A  CTL1  TO  0 
CONTROL  Gpio,2;0  (SET  CTLO  TO  1  A  CTL1  TO  0 
CONTROL  Gp i o , 2 2  (SET  CTLO  A  CTL1  TO  1 
(END  RESET 


870  RETURN 
876 

1180  Test  data:!  Generate  test  data 


•PRINTE&  IS  70S 
PR  IN’  C HR'S  .  i;  ' 

DISf 

BEEP 

INPUT  *What  is  tb»  Envslor  P»a*  A»r 1 i tud»?* .Epj 

BEET 

DISf  ’  GENERATING  DATA- 

P«ak_vaiue_da*a-r 
Ga l_cyc_f rg*200 
Ab*  0 

ALLOCATE  INTEGER  N»»<0:N«,*a  U 
FOR  X-0  TO  Ns* A  -  1  STEF  A 

Ac*Epa*SlN< ( PI • ( ♦  Ab  > > / < G*Ga i_cy c  <  rg > *2 .  7*S IN ( ?«P I *Ab*Gf /S* 
FOR  Y • X  TO  X  ♦  ? 

Nnmr  (  Y  :  *  I N  T  ( Ac  >  ♦  1  ?8 
•PRINT  Y.tW<Y> 

OUTPUT  ^Buffer  USING  ** . B* : H»» . T • 

NE  X T  T 
Ab* Ab* 1 
NEXT  X 
DISP 
RETURf- 
END 


viKy 


•A179  742  LASER  BALLISTIC  SENSOR  DEVELOPMENT!!!)  BOEING  AEROSPACE  2/2 
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UNCLASSIFIED 


F/G  9/2 


NL 


1 

f 

& 

Envelop_f22 

552 

Epa 

Err  rout 

280 

288 

1188 

1196 

FI 

306 

<-DEF 

312 

321 

396 

404 

.« 

F2 

306 

<-DEF 

326 

335 

459 

467 

Fb 

267 

268 

316 

330 

•j 

F i le_array 

Fi le_data 

Filt  bu  int 

,v 

Filt  fi 

275 

411 

418 

419 

431 

519 

S20 

527 

Filt  f 1 1b 

“lit  fllr 

K 

-D  u 

1 

4-*  4-» 

— <  — « 

-a  -e 

U-U. 

Jw 

Filt  f2 

275 

474 

478 

479 

493 

541 

542 

549 

Filt  f 22b 

A 

Filt  f22r 

Filt  f 2b 

Filt  f  2r 

*>. 

Gal  eye  f rg 

284 

288 

1192 

1196 

Gt 

127 

288 

316 

330 

1196 

.V 

Gpio 

119 

252 

253 

255 

259 

260 

866 

867 

868 

•jt 

Hpib 

120 

Si 

I 

722 

723 

724 

725 

728 

729 

730 

731 

If  1 

271 

275 

314 

315 

392 

396 

404 

430 

431 

516 

* 

537 

538 

* 

;* 

If  2 

273 

275 

328 

329 

455 

459 

467 

492 

493 

540 

M 

559 

560 

J 

384 

41  1 

418 

419 

421 

452 

474 

478 

479 

481 

i  1 

128 

245 

275 

360 

388 

431 

570 

580 

r«s 

k2 

129 

454 

v 

K3 

130 

<k 

M 

395 

396 

401 

403 

404 

409 

458 

459 

464 

466 

jSji 

467 

472 

S' 

i 

V* 

Max_x_axes 

Mmm 

Mn 

Mult  fun 

286 

290 

292 

1194 

1200 

1202 

M 

268 

314 

315 

328 

329 

360 

395 

403 

458 

466 

•’*3 

(inn 

Ns 

Num  samp 

240 

241 

242 

243 

706 

796 

1194 

1195 

.•S 

IhS 

Num  samples 

239 

240 

245 

248 

258 

286 

287 

823 

828 

844 

'V 

852 

Op 

314 

315 

321 

323 

328 

329 

335 

337 

354 

565 

566 

575 

576 

633 

639 

676 

685 

695 

696 

695 

'll 

702 

Peak  value  data 

283 

1191 

Peak  value  fi 

385 

418 

419 

peak_value_f 1 1b 

Peak  value  fltr 

Peak  value  fib 

•<i 

1 

Peak_value_f 1 r 
Peak_va lue_f 2 

Peak _value_f 22b 
Peak_value_f  22 r 

453 

478 

479 

Peak  value  f2b 

Peax_value_f 2r 
F'roduct__t  1 

396 

397 

404 

405 

V 
:*' 

V 
w 

•:' 

V 

91 

kJ 

V' 

■>  -  ^  i’  r  "  r,  r  «>  n  r*  r 

■  ■  ’  ■  ,  .\V,  V  >:  V  v 

wrr*  or*  t 

•\\  ->  *'*Vi 

v*. 

'V 

»V‘i> 

Product_f 2 

459 

460 

467 

0 

375 

377 

388 

Ob 

371 

377 

580 

Cic 

536 

537 

538 

576 

580 

Occb 

243 

576 

Occr 

243 

566 

Qj 

Opsb 

645 

668 

675 

Opsr 

597 

623 

630 

Qptb 

664 

666 

668 

Optr 

619 

621 

623 

Oq 

365 

700 

703 

Oqq 

716 

724 

760 

Qr 

370 

375 

570 

Osb 

190 

646 

649 

Qsr 

189 

598 

601 

Ot 

180 

<-DEF 

181 

623 

630 

632 

Quot lent 

360 

361 

R 

717 

738 

745 

Ra 

605 

610 

611 

656 

657 

659 

Range_stored 

Rmdr 

702 

703 

Rob 

573 

575 

Ror 

563 

565 

Rot 

245 

Rotb 

243 

575 

Rotr 

r 

243 

565 

O 

st 

126 

268 

288 

Ssp 

Start 

361 

364 

370 

Sum_;f  1 

393 

397 

405 

Sum_f 2 

456 

460 

468 

Sys_del_bl 

134 

371 

Sys  del  red 

131 

370 

Ta_bl 

245 

685 

762 

Ta_red 

245 

639 

760 

Tab 

646 

649 

668 

Tar 

598 

601 

623 

Tb 

Test_data 

760 

762 

764 

T  i  c_m_x 

Tot_angle_peak 

353 

1  t 

X 

287 

289 

295 

458 

459 

466 

796 

798 

800 

Xx 

372 

374 

394 

Y 

Z 

hi 

289 

290 

292 

797 

798 

799 

*■  String  Variables 

rta$ 

213 

214 

224 

NameS 

742 

747 

818 

Titles 

754 

*  I/O  Path  Names 
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468 

454  706 


558 

559 

560 

563 

566 

570 

573 

677 

678 

632 

633 

678 

671 

673 

675 

677 

626 

628 

630 

632 

633 

716 

762 

774 

668 

675 

677 

623 

630 

632 

182 

183 

184 

185 

186 

187 

188 

668 

675 

677 

759 

773 

655 

612 

614 

617 

623 

630 

653 

662 

668 

675 

196 

371 

411 

474 

675 

685 

630 

639 

388 

395 

396 

403 

404 

423 

ASA 

467 

483 

758 

760 

762 

764 

765 

802 

804 

806 

1195 

1199 

1205 

457 

562 

595 

690 

293 

1199 

1200 

1202 

1203 

800 

801 

802 

803 

804 

819 

823 

825 

843 

846 

1 


eBuf  ter 

248 

251 

258 

292 

8S0 

852 

855 

1202 

•File 

825 

828 

846 

848 

§Gpio 

191 

258 

■*  Line  Labels 

A_to_d_reset 

198 

883 

<-DEF 

Arc_t 

588 

<-DEF 

Buf i er_array 
Bw_label 

220 

228 

791 

<-DEF 

Check_stat 

Check_stat_r 

Comments 

1 1 

<-DEF 

Comp_angle 

204 

345 

<-DEF 

Con_f 1 

380 

<-DEF 

Con_f2 

449 

<-DEF 

Constants 

Disp_con 

Disp_data 

Disp_data_1 

123 

<-DEF 

Disp_data_2 

Disp_out 

117 

207 

708 

<-DEF 

Dma_transf  er 

250 

<-DEF 

Doppler_cor 

51 1 

<-DEF 

Env  fl 

425 

<-DEF 

Env_f 2 

485 

<-DEF 

Err_rout 

113 

781 

<-DEF 

File_array 

225 

835 

<-DEF 

F l le_data 

219 

812 

<-DEF 

F i 1 t_bw_i nt 

221 

227 

265 

<-DEF 

Index 

81 

<-DEF 

Initialize 

109 

<-DEF 

Main 

195 

<-DEF 

Mul  tj_f  un 

201 

301 

<-DEF 

Num_samp 

215 

235 

<-DEF 

Range_stored 

199 

210 

<-DEF 

Test 

277 

<-DEF 

Test_data 

217 

1180 

<-DEF 

Var  list 

139 

<-DEF 

* 


Line  Numbers 


226 

240  846 

246 

249 

363 

372  707 

390 

624 

526  523 

540 

546 

548  545 

562 
571 
576 
600 
626 

634  599 

635 
648 

679  647 

686 
691 
716 
729 
736 
763 
816 

818  819 

820 
827 
830 
858 
884 
895 
902 
907 
911 
951 
1254 
1274 

-*  SUB  Subprograms 
Endcat 

Unused  entries  *  3 


602 

650 
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